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Wf&Sf  or  VAPOB  VELOCITY  OW  fOM  C0WBBMUfl3Cir 

by 


(Mmltted  to  tl»  D^partsMnt  of  Moebamical  Engiiioering  csa  Jamaary  XI, 
X955  In  partial  fulfillment  of  the  requirtBaents  for  tte  AiigrMi  of 
Master  of  Science  iri  Mecbanicti.1  Engineering.) 

In  a  vei'tical  plate  or  tul>e  ecxidenaer,  ttus  liquid  eoisdeoMte  layvr 
(filtt)  flows  dowi-MBrd  \»:ider  the  influence  of  gravity  foreea.    laaiinar 
flow  generally  eatiats  on  the  v^V^  part  of  the  surface,  but  turbulftat 
flow  will  exist  ori  the  loyer  parts  if  the  eoadftaaatlon  rate  is  suf- 
ficiently high.     If  the  vapor  velocitits  are  anpfraeiahle^  the  vapcar 
exerts  a  aiiear  stress  on  the  liquid  •^rap<»-  interfaec,  vhieh  aay  affeet 
tte  heat  transfer  coefficients. 

The  ease  of  Imtixmr  flow,  vsadier  the  i»fluttnoe  of  this  impor  slasar 
stress »  has  be^  analysed  by  Kusselt  (1).    By  assuaine  the  "uc^ivex^^al 
velocity  distribution**  of  Prandtl-KiJcuradse,  Ssban  (10)  obtained  ana- 
lytical predictions  far  the  turbulent  film,  b«t  neglected  the  effect 
of  vapor  velocities.    Utoia  thesis  eacteiids  the  previous  wcrk  by  tlevelc^- 
i£ig  eqfuatii»2S  lAiich  include  the  vapor  shear  stress  effect  on  a  turbuX<.^K.^ 
flte.    Two  criteria  far  transition  are  XHroposed.    Aj^Iication  of  t^iese 
criteria  penait  the  present  results  to  be  coabisuad  with  those  of  Husselt 
to  sive  analytical  heat  trazssf«r  coefficients  for  the  case  of  laminar 
flow  on  the  \x&per  portion  of  a  plate  and  turbulent  flow  on  the  loww 
portion. 

The  resvklts  of  the  analysis  are  presented  in  the  foxis  of  diarasion- 
lass  plots,  FipsTM  7  to  12,  of  heat  transfer  coefficients  versus  li<3Uid 
Reynolds  Htnber  at  the  bottoa  of  the  plate  and  also  versus  plate  l«igUi. 
Dlmensionless  shear  stresses  and  liquid  Frandtl  mxabers  are  pansistnrs 
in  these  plots.    The  analytical  results  are  shown  to  agree  closely 
with  Colbum's  correlation  e<^iation  (9)  in  the  range  IL    •  2  to  3. 
However,  considerable  deviation  exists  in  the  very  hig!n  or  very  lot?  W^^ 
range.    The  Pinndtl  Ifushers  covered  by  this  analysis  extend  from  0.01 
to  10. 

An  ejqperiiaental  set-up  for  investigating  the  hydrodynamics  of  a 
falling  liquid  film  is  described.     PreliEinaa-y  results  indicate  that 
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the  (mAlytieal  equations  are  able  to  predict  film  thlr  t^. 

eurately  fcK*  low  EejnaoMs  Numbers.     Hovevcr,  at  hl^Sae;  jcibers, 

BOBS  deviation  vas  observed.     It  is  conclwied  that  the  Prandti'^flkuredse 
velocity  dlstrlbutlcsi  Is  generally  adequate  for  the  present  purpose 
of  predlctljsg  heat  transfer  coefficients,  but  that  further  study  is 
necessary  to  verify  this  golnt. 


A  procedure  utilising  a  hot  vire  saMMaseter  is  dttMsribed,  vhich 
appears  to  hold  proalse  la  dettersiinlng  transition  points  in  the  Ix^uid 
fila. 
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TlM  iiifliWQ^e  Off  a.  v»per  sbeariog  screes  cjeai  c«ad«aciz^  f  11a  itieax; 
transfer  cocfflelenta  hui  been  ©feown  by  several  investigators  (1,  2, 
5,  4)  to  be  qttite  » JgniTieant .     Ho'wever,  thia  viariAble  ia  often  neglecte^i 
In  tfae  interest  of  a  iWTWiamhle  tbeoreticttJ.  oEuBlyBis  vith  the  result 
tlMit  experimentally  cletermlned  values  of  beat  transfer  eoeff iciente 
have  often  differed  froa  predicted  v^Llues  by  several  fold.     Tbe  devia- 
tion has  been  maob  nore  noticeable  in  t±ie  case  cf  vez*tical  tube  conden- 
sers where  tbe  cotidensate  flovs  in  film  fora  on  tiae  inside  of  the  tudMW. 
Tbia  seems  raascauftble  ^ibietx  oa&  ecs^iders  tbe  fact  tbat  high  vapor  vel- 
ocities are  more  easily  attainable  in  this  case  than  in  the  case  of 
horizontal  tube  ecndeosers  or  vertical  c<»»ien0ers  with  the  condensate 
in  tJiS  shell. 

Altboogh  higher  film  eorsdenaate  coefficients  can  be  obtained  vith 
big^  vapor  velocities  it  aust  be  rerssiabered  tbat  for  vertical  tube  con- 
densers, tl^  price  (^  this  Increase  in  heat  transfer  ability  is  %&  added 
pressure  drop  in  the  condensing  section.     In  the  case  of  steem  condeQ» 
sers  it  is  also  to  be  noted  tlmt  the  cooling  vater  side  coefficient  is 
usually  controlling  and  little  can  be  cpaiaed  in  overall  heat  transfer 
coefficients  by  decreasing  the  condensate  resistance  to  heat  flov.     itov- 
ever,  many  industrial  appllcaticms  involving  the  coasd«ru»ation  of  organic 
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vapOTB  are  concerned  vlth  corvdensate  coefficients  &a   high  aa,  or  hi^- 
er  than  the  cooling  water  coefficients.  It  is  in  eaaee  such  as  these 
that  reductioais  in  size  and  vei^t  nay  be  possible  through  the  Judicious 
use  of  this  vapor  shetur  effect. 

A  condensate  film  flovliig  vertically  under  the  ixifluetice  of  grav- 
ity and  a  vieqpaf  shear  stress  lends  Itself  reaiarkably  veil  to  theoret- 
ical attack.  It  is  the  purpose  of  this  thesis »  therefore^  to  attaspt 
to  further  the  previous  aoalyses  by  derivirig  express icms  for  average 
heat  transfer  coefficients  which  include  both  Imsinftr  and  turbulent 
eonAensate  l^styers  under  the  influexu^e  of  a  vi^or  shear  stress.  To  sizs.- 
plify  the  i:>esulting  expressions,  aiid  at  the  saste  tiae  give  results  of 
universal  applicability,  the  case  of  a  flat  plate  is  considered.  Tia» 
equations  nay  also  be  applied  to  ti^bes  as  loog  as  the  thickness  of  the 
e<»3densa.te  film  is  small  eoepared  to  the  tube  radius,  and  as  Icaag  as 
the  effects  of  a  preesta*e  gradient  may  be  neglected. 
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In  tfe€  following  ojoalysis  tlie  case  of  a  flat  vertical  p:u,te.  Fig- 
ure I,  is  Q<xmid9rt&*    Ccmdezisatlon  io  of  a  pure^  saturated  vs^r  witb 
no  tioQcoridfixisibLe  gMes  piresent*    AtmKqpftions  Inelxide  eosuitJBit  \iall 
t-p.«tur.  a:«l  constant  fluid  properti*..     The  specific  heat,  c^,  den- 
«ity,  p,  end  th«ra&I  condiiOtiYity,  k,  may  be  detensiAdd  at  the  arith- 
metic average  fil»  ten^rature.    Ihs  viscosity,  yi,  l»3wever,  should  be 
detenained  at  the  special  avera^^e  fllsa  temperature  reeotamended  by  Col- 
btam  (5),  and  defined  by 

The  vnpor  shear  stress,   'tV,  is  coasldered  as  being  directly  doimvard, 
and  being  constant  along  the  length  of  the  plate. 

The  laciiaar  analysis  essentially  follows  that  of  f!tu»selt  (l),  fOntre 
the  f  lis  is  In  streasiline  action  and  is  ui^uder  the  iiif  liience  of  a  vapor 
shear  stress,  ^^,  as  sha%si  in  Figure  1. 
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By  asfiuE&ing  st«ady  flov  and  neglecting  fflGiaent\&^  ehangfts,  a  force  b&l« 
for  tlae  8i3flded  elesteat  sho^oi  glvMi« 


T  -  -I-  (p  -  P^)  K  -  *)  -   T^  -  0 


(2) 


This  gives  the  Ilni^ur  distribution  of  shasr  s'trsss  slacrwa  in  Figure  2. 
Since   f  »  '^  at  X  «  0,  the  -Mall  flheor  stress  is  giveti  by 


'^»  -  -f-  (P  *  P„)  ^«  *    "^^ 


&. 


(5) 


The  effect  of  the  racmeaatum  crossing  the  eloneai  faces  vitifuln  the 
liquid  is  quite  aoell.  'JDo  include  the  nonefftun  effect  at  the  liquid- 
vapor  interface  the  tcna  %Aiich  should  be  added  to  Equatlcm  (2)  is  (l/s^)v„  v.  <^^ 

o  z,x^ 

This  represents  the  z  direction  velocity  of  the  vapor  crossing  the  li» 
VOld-vapor  interface ,  N^toich  is  very  nearly  equaX  to  -^ac   liquid  surface 
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velocity,  v        *  and  not  th«  velocity  of  tl»o  r&ijp&r  flow  in  ti»  mala 

*  o 
wver  strMBi.    n&is  mcBMrntvci  taxm  «lfto  is  tMUAlly  qiaite  flgnall. 

FroB  ltevtoQ*s  dvfiniticn  of  vlscofllty^  we  }m.ro  tor  laminar  flour, 

Si&8tituting  in  Eq]uation  (2)  for     i    glTM, 

This  expr^MioA  HUiy  be  integrated  acz^ose  the  filia  vitli   LAr  eaoMtazit , 
giving  the  velocity  profile  as 

^  -  "f-  (P  -  <^v^  ^V  •  ^'^/^^  ^  "'^' "   ^  *  ^^> 

Dm  Bass  rate  of  f  lov  per  \aiit  width    oi  the  vertical  plate,  P   , 
■ay  be  found  by 

C&rryijag  out  this  integration  gives    I      in  terns  of  filas  thielmess,  x  , 


s(p  -  P,)    -|-  *    ^a     °    I  (8) 


Ifiaation  (S)  nay  be  differentiated  with  respect  to  x  to  give  a  relA* 
tiooship  between  the  change  in  mass  flov  rate  per  unit  width  and  ths 
change  in  filja  thickness, 
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g  (p  -  Py)  x/  .  ^ ^-        dx^        (9) 


A  linear  teaperatiare  gz^adient  laay  be  assuBied  to  exist  within  tha 
coQd<aisate  fila  without  introdiicing  serious  error.     Befenrence  (6)  gives 
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a  coriplete  tre«ittaent  of  tbs  SMi^itud«s  liivolved  in  this 
Coxiaidsr  the  coatrol  volume  shown  in  Figure  5. 
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If  the  enthalpy  of     '     entering  throijgh  the  \9per  surface  of  the  con- 
trol voXvBte  is  taken  at  a  s»en  film  tes^erature.  It  can  be  aeesa  tiiat 
tte    1      leaving  throu^  the  lower  surface  of  the  control  "s^ltKae  vill 
fldso  be  at  this  sane  f  iln  tes^r attire  and  hnnce  will  not  enter  isto 
«a  energy  balance.     It  is  also  evident  that  the  enthalpy  chSB^i  tst 
the  incretaent  of  condensate,     J '     ,  will  be  equal  to  the  enthalpy  dif- 
ference bet^men  saturated  vapor  and  a  liq^uid  at  the  mean  film  tempera- 
ture.     I^ie  nmy  be  designated  h^    and  has  been  evaluated  by  Colbum, 
Millar,  and  'Westwitcr  (?)  as  hi    «•  h^    *  5/^  c    ^t.     On  the  basis  of 
the  previous  discussion  SA  energy  balance  for  the  control  <^olvae  shoua 
in  Figure  >  sLsy  be  aade  by  equating  the  heat  transfer  to  the  change 
in  enthalpy  ot  the  inereraental  increase  in  condensate,  hi    iC , 
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-^  dz  -  h^^  cj  r  (10) 

Since  the  fili^  considered  ia  in  Laaiinar  aotioa,  the  heat  trans- 
fer coefficient  isa^  be  talsen  as  the  ttoertaal  conductivity  of  the  fluid 
divided  by  the  filn  thlckneee, 

-i-  .  h  At  »  —  ^t  (11) 

A  x^ 

Eliminating  q/k  betve«i  Eq^iations  (10)  aod  (11)  tmA  aubetituting  for 
i.  ^      its  equivalent  from  Equ&tloa  (9),  ve  have  a  differaitlal  exprea- 
ftioQ  ooQtainin^  only  dz  and  dx    as  variables^ 


kAt     ., 

-r-l dZ 


^  p(p  -  p^)x^^  *  ^  p  r.x/ J  to^     (12) 


Equatic©  (12)  may  be  Integrated  frcm  the  top  of  ttoe  plate  vhere  x  «  0 


vhen  2  »  0,  to  a  point  x  ,2.. 

^  5 

"bj —  2  -  -a-  p(p  -  p^)  -np  >  -^  0  L^  ~Y-  (13) 

By  reaorraaging  and  8i;^>stitutiag  t^  fraMtl  miaber,  N      <•  c  iVl^»  ^b* 
resulting  expirees icm  is  obtained. 


N.    ^.i^^^^^j^i^^'     ^  ^ ^^-"^^  (;f7(i.) 

The  following  diimmsialLess  lengths  may  be  defined  as 
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Further  dsflrie  a  dlmensioiileas  vapor  shear  stress,  L^  ,   as 


3rAAJg9 


(17) 

Substituting  these  definitions  Iri  Equation  (I^)  gives  a  dla^Ie  arid 
eonvenient  escpression  for  the  variation  of  filsa  thickness  vith  leu^jth 
of  the  plate,  z   . 


»*  .  (x/)^  +  V5     "T-.   (x^*)5  (18) 

Tbe  iu»efulnes8  of  the  dlBensionleas  quantities  defined  by  Eqaa» 
tlons  (15),  (1^),  an<i  (17)  aay  be  more  fully  realized  by  deriviiig  ex- 
pressiOiis  for  the  film  Beyxiold's  ^utaber  auod  meati  heat  transfer  coef- 
ficient aa   follows. 

Define  a  film  Beynold*s  RmbtKr^  !L  ,  by 

By  making  use  of  Eq:uation  (8)  for  P  ,  ve  may  write 

Bearranging  the  above  expression,  ve  have, 

^  ^  *  (21) 

Kov  by  substituting  the  previously  defined  values  of  x   and  C^ 
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l^e  origiiml  BtiwMlt  relation  gave  only  point  ^mXas6  ot  the  heat 
trmnafer  coefficient  and  mbb  difficult  to  «i{>ply  in  «&  actual  deolgn 
stuly  «e  it  invxjived  a  functioa  of  the  ecmdenaate  flow  rate*    Thle 
difficulty  can  be  avoided,  however,  if  a  mean  or  average  l»at  trans- 
fer coeff iclaat  is  defined  as, 

h     s       2 — TS,..  (23) 

By  stabstitutlng  Equation  (8)  for    P    in  E<iuation  (23)  «"8  heem 

{2k) 
Introdi2cing  the  Proiidtl  HOibar  aad  raorrangioi^, 

(25) 
Mftking  use  of  the  diffienaiooleas  tenos  z  ,  x    ,  and    t><r  ,  this  reduces 


(26) 
Eqiaationa  (18),  (22),  and  (26)  coetp^iftely  define  the  reaulto  for 


the  laminar  region*    For  given  values  of  the  dlDensionlese  shear  aferaaa, 
l^    ,  either  the  mean  heat  traiasfer  coefficient  group,  or  the  H^ 
can  he  plotted  as  a  function  of  the  dimenaloaaleaa  length  of  the  plate, 
%  *    Thie.  Involves  assuitlng  values  of  x      and  than  solving  for  s    trcm 
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tlon  (18),  M,^^  from  Equatlcai  (22),  and  kn  f^)     fron  JBqaatloa 
(S6).     This  has  teca  dontf  for  v«rlou»  li      aofX  farms  t^  b««ls  for  tlw 
curves  to  the  left  of  tlje  traasitiou  poiato  in  Figures  7  to  12.     The 
WM  of  these  curves  in  a  ilesi^gi  ease  vill  be  discmssed  la  a  Later  sec- 
tion, vfcere  ^»  results  of  &  tm^mOimt  f  Ito  as  wU  as  the  Imicar 
fila  ham  be«a  deriwd,  mA  are  plottM  in  tewue  ©f  ttese  sMse  ^ari- 
s^Xes. 


-.1  '5_ 


simoK  c  *  TRMsvnm 

For  A  liquid  floving  in  film  fonB  40m  «  Tertleal  •urfaee  it  9tmm 
iTOTable  to  postulate  that  there  vill  be  transition  froa  ligainar  flov 
to  turbulent  flow  at  a  critical  VeynaoMB  nvmbmr,  in  ataeh  the  sane  »«&- 
!a«r  that  transition  oceors  in  a  pipe  full  of  fluid.     .«1ien  \m  speak  of  a 
critical  Heyooldi  nvnber  it  is  to  be  renmatered  that  this  la  not  si  def- 
inite ^  distinct  point  at  >«tdoh  transition  vill  occur  in  all  eases. 
Rather  it  repres«tits  sane  sort  of  a  aean  value  vtoere  the  turbulfficuite  is 
sufficiently  developed  to  fflaterially  affect  the  velocity  profile  within 
the  f ilffl  and  vhei^  turb\ilent  velocity  fluctuations  will  be  evident  to 
tte  extent  at  affecting  sitch  characteristics  of  tim  tih&  as  the  heat 
transfer  coefficient. 

Before  proceeding  vlth  a  discuss  ic«i  of  trsuisiticm,  it  becooi^t  ne- 
ewaary  to  define  a  Reynolds  number  to  give  a  ccnaon  basis  of  cospari- 
son.     In  the  case  under  consideration  it  has  beeoBie  accepted  to  define 
tte  Reynolds  nvaSaer  in  terms  of  the  t&ass  rate  of  flow  ^r  unit  cross- 
sectional  area,  (v/S),  the  avera^  fluid  viscosity,  ti,  and  an  equiva- 
lent hydraulic  diameter  caeaxml^  tajcen  as  (ks/Y),  vliere  Y  1a  the  unit 
vldth  of  the  plate  dovn  vhioh  the  fluid  is  flowing.     Xn  tenat  of  these 
fuantities,  the  Reynolds  taxabmr  oay  be  defined  as, 

^  .  c:^ ,  CM!!i) .  jr 

Although  there  has  been  little  experiii^sntal  vork  vltia  the  specific 
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alm  of  determining  the  transition  point  of  a  fre«ly  falliag  vertical 
fila^  certain  points  zaay  be  inferred  fraa  taeat  transfer  data.  In  ref- 
erences (12)  and  (l^)  on  the  basis  of  sonevtrnt  limited  c'dta  it  is  rec* 
ontnded  that  a  Reynolds  nunber  of  ISOQ  ha   taken  as  the  transition  point 
for  the  case  of  aero  vapor  shear  stress.  It  vill  be  shovn  later  that 
considerable  doid)t  exists  as  to  this  vttl\»,  however  froa  the  present 
extent  of  knowledge  no  taare   precise  valiie  can  be  predicted.  Therefore, 
for  the  p\irpoec8  of  the  analysis,  this  criterion  for  transition  vill 
be  adopted.  One  of  the  objectives  of  the  eacperlmKital  work  undertaken 
as  part  of  this  thesis  will  be  to  determine  more  adequate  inforsaation 
on  transition. 

Although  we  have  adopted  a  critical  Reynolds  number  of  IBOO  for 
sttro  vapor  shear  stress,  there  reaaaliu  the  ccmslderi^le  problem  of  con- 
sidering how  a  shear  stress  oay  nodify  this  value.  Carpenter  (3)  has 
iRigcested  tiimX  under  the  influence  of  a  vapor  shecur  stress,  transition 
will  occur  at  a  given  value  of  wall  shear  stress.  ThXe  may  be  sjore 
readily  visualised  by  considering  the  layer  of  fluid  next  to  the  wall. 
In  lamlntir  f lov*  the  sole  effect  of  the  outer  layers  <xi  this  inner  layer 
will  be  a  shear  stress  exerted  cwi  it.  In  effect,  all  the  inner  layer 
''sees"  of  the  outer  layer  is  a  transmitted  shear  stress.  If  we  first 
eonsidsr  the  case  of  zero  vapor  flow,  but  increasing  liquid  flow,  this 
iansr  layer  wsxt  to  the  wall  will  trip  over  to  turbulent  motion  vknea 
the  Reynolds  nisaber  attains  the  value  1600.  ^S/cmemr,   all  the  imier 
layer  "knows"  of  this  condition  is  a  certain  value  of  tranai&itted  shear 
stress. 
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RbVjCOiMlderlng  the  case  of  a  finite  v«,por  shear  otress,  it  is 
rcMiccMble  to  as«me  that  the  layer  oeiet  to  the  mXl  vill  change  to 
turbulent  flov  \ihttnieyer   this  "critical"  ehoar  stress  is  exerted  on  it, 
vhether  by  o\ater  layers  of  li.-^uid  or  by  aa  outer  layer  of  vapor.  Tak- 
lag  the  wall  shear  stress.  To  ,   as  the  iac«t  oonvealeat  to  which  to  ap- 
ply this  critericai  of  critical  shear  stress,  ve  may  find  its  ^alue  fron 
the  known  cass  of  transition  at  a  Keynolds  number  of  iBOO  with  zero 
va^por  shear  stress.  Frca  Equaticai  (5),  T'ocan  be  expressed  in  tssms 
of  x^  for  r^  .  0, 

o 
Similarly  P     ,   snd  hence  Reynolds  naaSaer,  hP    ,   nay  also  be  expressed 

in  tenu  of  x_^  fro.  E»«tloo  (8), 

(^)tr«  -  1B<»  -  V5  -^  (  1  -  P>  )  x^' 
Eliainlkting  x    between  the»e  tvo  equations  ve  are  able  to  arrive  at  an 
•aqpression  for  the  critical  wall  shear  stress,  designated  by  (   T*©   )+«_^» 

(  ^'  >tr«  -     37  (^-  ^--^  L    (1  -  '''//'     J  (28) 

For  the  aore  gsxieral  case  where  T/*-  jt  o,  we  ai>ply  our  postui-itioai 
that  traasiticsi  vill  still  occur  at  this  value  of  a  critical  wall  shear 
stress  derived  above.  Prom  e^i:uation  (5)  we  have  an  expression  fear  l^o 
in  the  presence  of  a  vapor  8h»ur  stress.  By  introducing  the  viiue  of 
wall  shear  stress  obtained  in  Squat icm  (23)  and  rearranging,  thee  follow* 
\x%  eicpression  for  the  dlxaensicmlsss  film  thickness  at  transition  can 
bs  obtained. 
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Mbfltlt\xting  this  mlue  la  Equation  (£2)  giws  «a  eqimtlcai  for  tiw 
tsnuMtlticn  RcynoMA  nwrtbwg  In  tarms  of  thA  oasdy  variable,    Tl^    . 

^   ^     >traa  "  ^®^  *  2^(XV^)^/^T/  t  0.667(1  -  p^p)  (  '^J')' 

(30) 
Thle  funetlOQ  i»  plotted  in  F:Ue>»^  ^  ft^'ow  »  value  of  TV  «  v*)  to  ^^   «  9.2, 
Inspection  of  Bquetloa  (29)  will  reve&l  tbat  ^^heitevtii  X* >  -=^^ — r^-^, 

(iV'p^ 

tbe  diiaeiialoaleas  film  tbickr*«8»,  x     ,  vlli  be  nepitiy«  .at  trains  it  Ion . 
CSbviouely  thle  i«  lapoeeible  phyeic&lly  and  repreeente  a  li«alt  to  the 
applicability  of  the  previo\ysily  diecuseed  criterion  for  traaaition. 
In  the  hlglier  rent^see  of  vapor  ebear  a-brese,  v«  suiet,  tberefoz*e,  develop 
a  aev  ei'lterlon. 

For  the  hli^r  values  of    '^a^  it  is  9;uite  poeeible  thftt  the  fila 
thlclmeee  at  trahsltlon  approaehes  a  liaitizig  value,  perhe^  that  of 
tlie  leMlmir  sublaym-  thielcoeee  which  exists  in  ordinary  pipe  fl&t/* 
Carpenter  (^)  mgeests  the  following  tssfinitude  for  this  thiekoessx 


CMlztaril^  a  Msgnltmflle  of  5  esq^resses  the  Isminar  svflblaoi^r  thielosws . 
Aetvally  sa  appropriate  value  1^  soneKlMHra  between  3  and  perhaps  10. 
The  value  of  ':.  is  arbitrarily  chosen. 

Xf  Sqpation  (31)  is  solv^  for  To  and  this  expresslosi  siAmtituted 
ia  SQMitioa  (3)>  the  following  cubic  equation  relating  fll^  thicJtsness 
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&Qd  impor  «be«r  stress  is  obtaxjaed. 

For  a  given  value  of    'i^  j,  x      »By  be  fcnaid  frcBi  this  ©quation  througli 
a  "terifll  and  error  px-octsBB,    WitJi  this  valxie  of  x      aod  the  gives  value 
of     t  *-    ,  the  Beynolda  nistber  may  be  tovead  through  the  use  of  Fquation 
(22).     The  resulting  plot  is  sbowa  in  Figure  k  for  vsilues  of  'T--  >9.2. 
It  cafi  be  seen  that  the  assuEspticxi  of  transition  at  x      «  6  is  alioost 
equivalent  to  the  assynption  of  trieu^aition  at  a  Reynolds  mstber  of  about 
70  to  8o,  there  being  little  variation  erf*  BeynoldB  nutaber  with  vapor 
•hear  stress. 

This  appears  on  f  Ixiit  ob:^iervatlon  to  be  a  somewhat  startlngly  low 
Reynolds  mKbcr  at  \i^ich  to  expect  transition,  hovever,  there  is  a  log- 
leal  basis  for  its  selection.     There  is  very  little  good  data  availabler 
to  establish  conditions  under  whieh  transition  occurs  vhen  there  is  a 
si^lf leant  vapor  shear  stress,    A  single  point  can  be  obtained  from 
the  enthanol  data  of  Carpenter  (3)  by  plotting  his  experimental  points 
and  noting  the  regicm  of  departure  from  the  Htasselt  laminar  prediction. 
This  transition  point  at  a    '^   *»  110  and  Beynolds  n\saber  of  96.^  agree© 
fairly  well  '**ith  the  trans iticm  Beynolds  nusiber  of  71 '8  predicted  ujjder 
the  present  analysis.     In  tersas  of  x      the  transition  data  have  the 
value  of  6.9  which  is  close  to  the  6  of  Equaticaa  (31). 

Although  the  criteria  for  tr^nsiticai  as  portnayed  by  Figure  h 
have  not  been  adequately  sxibstiintiated,  the  results  are  nevertheless 
quite  reasooable  ood  are  Fart>bably  not  seriously  in  error.     This  entire 
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phsnaatonaa  of  tranaltloa  aa  affected  by  vapor  shear  atress  inust  b«  In- 

Testlgated  by  obtaining  much  more  experlrawital  dat«,.     In  epit«  of  lack 

of  adeaua-U  »ubstantiaticxi  it  8«eai8  reasonable  to  select  tke  traasition 

I    I'- 
crit«ri*  aa  ahova  la  Figure  '*  for  us«  in  tbe  ijresent  ai^^^^* 
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SifeWftS  D  »  TUPBULEHT  ANALYSIS 

When  condetisate  flov  rates  becoae  Large,  it  se«u)  reascmable  to 
predict  that  turbiileat  flov  viU  restilt,  with  considerable  effect  en 
the  heat  transfer  coefficients.  It  was  observed  experi»eatally  by 
Kirkbride  (2)  that.  In  the  ease  of  vapor  condensing  on  the  outside  of 
vertical  tubes,  laeasttred  heat  transfer  coefficients  differed  consider- 
ably from  those  predicted  by  the  Nusselt  laminar  analysis*  In  theae 
investigations  there  was  evidently  no  i^>preciable  va{K>r  stoar  stress* 

Colbum  (6)  reviewed  the  results  of  Kirldoride  and,  by  assuaing 
that  the  sole  reaiataaee  to  heat  transfer  lay  in  the  leainar  sublayer, 
was  able  to  derive  heat  transfer  equations  for  a  turbulent  condensate 
Isyer*  The  Prandtl-Nikuradse  velocity  distribution  vas  assuBftd,  and 
by  selectir^g  the  thickness  of  the  laminar  sublayer  at  about  x  «  10, 
a  satisfactory  correlation  of  the  data  i«s  obtained*  This  equation, 
as  presented  by  Colbiim  (8)  is 

\^  m  -  "^' ("-^  (f) 


It  is  to  be  noted  that  the  Prandtl  nuoBber  enters  here,  ^ile  it  did 
not  in  the  laminar  analysis*  This  does  not  seem  imusual  when  one  con- 
siders tlike  correlation  equations  for  beat  transfer  in  pipes  with  no 
phase  chao^ .  However,  it  should  also  be  noted  that  this  correlation 
was  applied  principally  to  Kirkbride 's  data>  with  only  saall  variation 
of  N_.  in  the  range  of  2  to  5* 


-^'±- 


S^BU^SiLl3&SL  i§ 


no  ld»tt«  9jLdart9btmao9  tilM  ^iimmrt  lLl%t  volt  iamlu&iMS  SmdS  S9tb&rtq 
i  uo  mL^  no  ai!ttiTftnfjT»  -io^v  to  MAO  •<(#  at  ^^mdS  (S)  •Jkl:ccWiUI 

.—A-h^u  -iMsuii  ioqj0r  mltkUtmqn^  on  v^«^£aii4»iv«  anr  iitrit  aooUiq^J^aimii 
aoUvMM  x^  «fta«  •MtditxlX  to  •ilmtri  mU  bmemirwrg  (6)  erwdloO 

»i>— mUbmo>  tamlu&wi  m  tD\  %aotfmip0  ttitmamtt  trntd  arlvtib  o^  mlifB  mv 
ItfHi  ,6s«nw«  Mw  cioilWl'x^alib  \  v  MtevnliM-IllMWVl  «IT    •ia%al 


(li\   Uv\)  ^-^-.^  =.  (^A\-cL 


MJ>  IJt  sXiiir  «Mflwi  n*txMi  tb^wr  XttovA  «tf  tJUlit  bm^cn  mT  ot  U  #1 

on  dttv  Mtqi  ilttors^  ^««(  -sot  ^aolSmuQ^  aolfml^'^i'/^"-  9df  ars»bk% 

s»l;rAlimoo  SilAi  Jmd^  Jb*tos  W  omSM  tduoiim  St  ^tmnmoiL    #«n>ai*n 


to  mi^iun  9dt  at    .N  %o 


fte  above  analyals  vae  extended  by  Carpenter  aod  Colbum  (3>9) 
to  include  the  effect  of  a  vapor  shear  str«8S*  In  i*«ference  (j5)  beat 
transfer  coefficients  vere  in ves tinted  in  a  5/Q  iach  copper  tube  9 
feet  long.  Special  effort  vas  aade  to  accurately  lasamire  pressure  drop, 
fr-om  vhich  the  magnitude  of  the  vi^por  shear  stress  could  be  obtained* 
As  the  result  of  the  experimexital  vork  aod  theoretical  investi^tions, 
the  folloving  correlation  equation  for  point  values  of  heat  transfer 
coefficients  vas  obtained: 

JUL    =  o.o'fs  [Nprj       T^ 

By  integrating  over  the  length  of  the  tube,  an  equation  for  mean  beat 
transfer  coefficients  was  also  obtained 

By  the  introduction  c^  the  dlu^nsioalesa  vapor  shear  stress  defined 
in  Equation  (17}>  this  equation  may  readily  be  converted  into  the  form 


k^ly^y^  o.o^^  (N,y  (r:)' 


(56) 


The  question  of  properly  evaluating  the  v«4>ox  shear  stress  beccoAS 
somevhat  complex  vd^en  it  is  realized  that,  as  condensation  is  taking 
place  along  the  length  of  the  tube,  the  vapor  velocity  is  continually 
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deereasing.     For  convenience  in  making  calculations.  Carpenter  and  Col- 
burn  reccsosaend  tsAing  a  mean  vapor  shear  »treas,  evaluated  from  the 
usual  plots  of  friction  factor  for  caae  ccxjqponent  f lov  at  a  va{>or  atass 

velocity, 

/     2  2  \  ^/'^ 

^"  I 1 7 

(:57) 

where  G.  i»  the   value  of  G  at  the  top  of  the  tiitoe,  and  Gg  is  its  valu« 
at  the  bottom.  It  is  felt  that  this  represent a  soaevhat  of  an  overs ia- 
lificatlon  and  that  the  friction  factor  should  at  least  be  evaluated 
froB  the  limited  data  cm  cc»icurrent  tvo  phase  flow  available  (11).  In 
reference  (11)  it  is  showa  that  there  is  coaisiderable  variati<Mci  of 
friction  factor  from  the  smooth  pipe  data,  mainly  due  to  the  fact  that 
the  liquid  layer  presents  a  rough  surface  to  the  flowing  vaipor. 

One  point  especially  worthy  of  coraowit  in  connection  with  the  w>rk 
of  Carpenter  and  Colbum  is  their  ec»isiderati(»i  of  the  aoBimatua  fltix 
across  the  liquid -vapor  interface.  In  their  derivation  they  evaluate 
this  term  as  the  product  of  the  ir»:rement  of  condensate,  <i '     ,   and  the 
main  stream  vapor  velocity.  This  in  effect  assuaes  one  dieiensional 
flow  of  the  ra^oT,     In  view  of  the  fact  that  a  vapor  shear  stress  is 
also  considered,  it  seems  scmeidiat  contradictory  to  assume,  ^^ben  eval- 
uating the  mcnentun  tens,  that  a  zero  velocity  gradient  exists  at  the 
interface,  and,  when  ecKisiderlng  the  vapor  shear  stress,  to  assune  a 
finite  velocity  gradient.  It  SMus  more  realistic  to  visualise  the 
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layer  of  viqpor  iiaaediatcly  next  to  the  interface  as  flowing  -with  essen- 
tially the  Moe  wlocity  as  the  interface  itself.  In  this  caee  it  can 
be  reeuiily  seen  that  such  mc2iaentuEn.  terms  vill  be  aaall.  It  is  through 
this  line  cxf  reasoning  that  this  mociwitua  tera  bas  been  neglected  in 
the  present  analysis. 

Quite  recently  Seban  (10)  has  applied  analogy  caJ-cui^tiona  of  the 
Fraodtl-Karmtm  type  to  a  turbiilent  condensate  layer,  assizaing  the 
**  universal  velocity  distribution"  of  Frandtl-Kikuradse .  The   improve - 
sent  over  previoxAS  analyses  lies  in  the  fact  that  the  resistance  to 
heat  tr«uasfer  of  the  ccmplete  coodensate  layer  is  considered,  i*ather 
than  jUst  the  resistaiuse  of  the  laminar*  siablayer*  Seban*s  analytical 
results  extend  over  the  considerable  ran®e  of  H-^  of  0.01  to  5  and 
a^sree  veil  vith  the  asipirical  results  of  Colbum  (S)  in  the  raaga 
Hp^  «  2  to  5«  In  the  derivation  Beban  ne&Lected  the  jK>ssible  effects 
of  a  vapor  shear  stress.  It  is  the  purpc^e  of  the  present  analysis 
to  extend  £cban*s  results  to  include  this  effect  and  to  cosplete  the 
story  by  obtaining  mean  heat  transfer  coefficients  which  include  a 
ls«1nar  region  at  the  top  of  the  plate,  followed  by  &  turbulent  rarige 
over  the  lo\mr   regions. 

Turbulence  has  no  effect  cai  the  force  balance  for  a  control  vol- 
VB»   in  the  ccsidensate  film  if  ve  think  of  the  shear  stress  as  an  appar- 
ent X   rather  than  the   KevUaiian  viscous  shear  stress.  Therefore, 
the  resulting  force  balance  equation  will  be  identical  with  that  for 
the  lazainar  ease,  and  is  given  by  Eqpiaticais  (2)  and  (3).  The  saoie  lin- 
ear distribution  cf  shear  stress  vithin  the  film  vili  exist. 
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To  -   -^  (P  -  Pv'  *o   *    ''^•^  <') 

3ine«  Nevto»*s  Isv  cdT  viscous  shear  vill  not  hold  In  the  tuEl^ulsnt 
regiod,  it  is  HapcAsible  to  obtain  aoalytieally  sn  «s|KrwM»iai  for  tiis 
Telocity  distribution  vithin  the  coadsmuite  film.    This  szialyais  viU 
be  baiwd  Goi  the  assumption  that  the  velocity  dlstrlb^ion  vill  be  that 
of  the  FraMtl-MiiBiraaM  formulation  as  given  in  refercmee  (12)  . 

For        0   <  X*   <   5  V*  •  X* 

5    <  x"^   <   30  v"  •  -  5.05  ^  ;2.u  in  x"^ 

50  <  x"^  v''  «  5.5  ♦  2,5  In  x"^ 


(53) 

The  apparent  shear  stress,    '^    ,  for  turbulent  flow  m$af  Imi  esprsttseA 

as  ttie  SUB  of  the  viscous  shear  stress  S£td  U^  edsy  dtiffufliirlt^  i>'me^ 
the  velocity  gredient. 

tfiiere     6k*t  is  the  eddy  diffusivity. 
T^carraagias,     ,?9  X  -    fy^  +  £^^    ^ 
Mbstituting  for    T     frctt  Eqiuatioc  (2) 
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ConBlderinfi  the  region  in  tha  film  'idsere  30  <   x*<  x^*",   it  may  be 
aasiaaed  tiiat  lJ<<<Ey^,    Solving  for   <£»*>y^    , 


6^    - 

F 

Jx 

(^) 


Itje  term  dv/dx  nacf  b©  evaluated  for  tiae  region  under  coosideratioo 
(  30  <  X*    <    X  * )  froa  the  asBUfted  universal  velocity  distribution, 

v**"  «  5.5  4  2.5  In  X* 


dv/dx 


2.5 
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IT      I 


/         X  (i^l) 

tl»  contrcl  vt>l«fic  for  the  energy  balance  is  showa  in  Figure  5» 
Also  shown  is  the  type  of  teaperature  distribution  to  be  eacpected  for 
tiai^ulent  flov. 
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Por  the  s.T««ied  coiitrol  volume  it  can  be  »e«r.  timt  th©r«  vili  be  iio  ea- 
tnalpy  cnangc  for  trie  portioa  of  coixdcasate  entering  trirougii  tae  upper 
surface  asd  leariiig  through  tae  lo^wer  surface.     Further,  it  is  appareat 
ti-at  tae  eiitiiaJLpy  eJatttige  of  the  increment  of  coodeiisate  will  be  esseri- 
ti&lly  ttJAt  from  aatxarated  vapor  to  satvurated  liquid.    Therefore,  «ti 
eriergy  balance  will  be  of  tirf»  f ora, 

-9L  dz  -  h-      d  r 
A  fg 

Aa  tiie  left  side  of  tbe  control  surface  is  moved  across  t^  fiJjl 
tovard  the  wall,  the  above  statemeats  viU  hold  true  except  for  the 
thia  layer  ««rt  to  the  -wall  where  the  teaiperatux-e  distributioa  experi- 
ence* a  sharp  ctaaai^  in  slope.  Since  this  laorer  is  thin,  tae  effect  of 
sia>cooltag  Bwy  be  neglected  for  the  turbulent  f  11»  and  q/A  will  equal 
(q/A)  ,  c^  the  heat  transfer  per  unit  area  at  the  nail.  Introduciag 
this  into  the  energy  balatice,  ve  have 

^  A  ^o    A      fg  dz  ^   ' 

The  equation  for  heat  transfer  across  a  turbulent  f ila  may  be 
written  in  the  fois 

i^/Al  =  (a  +  £.)  i± 

u 
vhere  a  «    -^      «  therraal  diffusioa  coefficient 

£;^  «  eddy  diffusion  coefficient 

Through  the  application  of  the  heat  •  aoBentum  trar^sfer  analogy  that 

6l^  »     6w,     ,   aftd  knovlng  the  velocity  gradients  throughout  the  film 
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frcB  tbe  asstned  Telocity  distribution.,  ve  are  able  to  obtain  an  expres- 
8ioa  for  the  teaparature  gradients.  By  substituting  <fi)  for  St»,  in 
Bquatiou  (40)  and  tbe  velocity  gradient  from  Eqiiation  i^l),   ve  obtain 
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Substituting  this  expression  in  Equation  i^^),   the  desired  eiqpr^ision 
Involvlog  dt/dx  and  x  as  the  only  variables  is  obtained* 
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Solving  for  dt  ar«d  introducing,  for  coovenience  salte,  the  foILoving 
notatio»,  J?  =  —  =  J  X  To  ^  ^  V^  ^  ^fi'-^/zJ^ 
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Tias  condensate  fila  may  be  divided  into  the  folloviri^  regicme: 
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X,   -^ 


X'^=S'       x+  =  jo      x*  =  yo*' 


FIGUEE  6 

Equation  (H)  contalna  oiiXy  tvo  variabl««,  t  axid  x,  since  ve  have  shovn 

tkiat  q/A  will  be  constant  across  the  film.  This  ttquation  may  then  b« 

integrated  between  the  lis&its  t           to     t  and  x,   to  x  .     After  con- 

X,  X_  X  o 

X  o 

siderable  al^braic  saaxiipulation  atid  siaqplification,  tine  re   :!lt  of  this 


integraticm  across  the  filA  aaay  be  expressed  as. 
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Equation  (^5)  gives  «a  express  ion  for  the  t««Bp«erature  droj)  ocrosa  the 
highly  turbulent  core  of  the  condeosate  fila.  By  utiliricg  the  results 
of  a  previous  analysis  by  Von  Kanooji  (l4)  foi'  the  temperature  difference 
between  the  vail  and  the  outer  edge  of  the  buffer  layer  (x  »  50),  we 
caii  obtain  an  expression  for  the  total  temperature  difference  frota  the 
liquid -vapor  interface  to  the  vail.  Vou  Kamaci's  results  aay  be  put 
in  the  form 


(k6) 


Addiiig  Equations  (43)  t^id  (^)  to  obtaisi  the  total  teraperature  drop 
across  the  f lln. 
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Let  the  right  hand  side  of  Equation  (47)  equal  Fg,  thea  this  ex- 
pression may  be  written 
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For  X       >    50 

o    ^ 


p  =  /^'(-o^  -»-  3k^*  -f  2.$x^*  In  x^*) 


For  5  ^  \<  50 


r  =  /;/  (12.55  -  8.05  X  *  V  5  X  *  In  x  '^) 


Different iat lag  with  respect  to  x      gives. 


For  x^     >    30 

o     ^ 


^^,  '/-^(S.S-a.Slnx^*) 


For    5  <  x„*  <    50 
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For  ccaavienience  define  F-  aa  follovs. 
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(^) 


Jy/     ' "  '  -^^^  -  -'  *"  -o  '  (50) 


For  x^  >  50        F^  s   5.5  +  2.5  In  x^^ 
For  5  ^x^''^^  50    Fj^  =  -5.05  +  5  In  x^""  (51) 

Tben  Squat  ion  (50)  laay  be  written  •• 
For  X  *  >  30       -4^^  =    f^^i 

For  5  <  x^*  <  30   -^^   =  /A^  i'T 

(50a) 
Eguatioa  k2  gives  the  results  of  tlis  eiMergy  balance  ia  the  form 


^  A  'o        "fg    j2 


This  may  be  written  as 
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fktai  Equation  (^7a}  ve  have  an  express Ion  for  (q/A)  as  &  fuoction 


of  film  thickness,  vhich  m&f  be  rearranged  in  the  form 


f  _    At  /"  Cf,  U:) 


i\-^-/f. 


ni^ 


M  J 


(1^  A^ 


A 


/: 


(55) 


Since  it  has  been  s2u>vn  that  (q/A)  -  (q/A)^,  this  term  way  be  eliminated 
between  Equations  (52)  and  (55).  Substituting  for   7-+  froa  Equation 
(50a),  we  obtain  an  expression  containing,  as  the  only  variables,  x  , 


dx  ,  BxA  dz. 

o  * 


F,  r. 


(Xt) 


ti'* 


V/^ 


M 


(1-A//J 


Jx/  =  ^t 


(y(f)'^^ 


This  nay  be  integrated  bet>*een  the  limits  x      to  x      and  z^  to  z,  idi«re 

y^^  and  z^  indicate  the  Talues  at  the  traasiticxi  froai  lawinar  to  turbu- 
ot  t 

Lent  flov.    llws  resulting  sfiiatioii  «ay  be  put  into  the  convenient  fora 


(Z-Zt)  At   =■ 


n'/s 


M 


ui-m 


dx.-" 


c. 


V5 


^ 


t» 


/^"; 


(5^) 


Squation  (5^)  now  gives  the  desired  expression  for  the  distance  down 
the  plate,  z,  as  a  function  of  film  thickness.  In  solving  for  this 
distance,  the  value  of  x   at  the  transition  point  must  first  be  solved 
for  froa  either  Equation  (29)  or  (32)  and  then  the  distance  to  the  trans* 
ition  point,  2^,  found  from  Equation  (18).  ITie  distance  fr<aa  the 
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tranaitictti  point  to  the  point  under  cojuiider»tion  may  then  be  calculated 
frooi  x^iquation  (5^)  by  evaluating  tiae  Integral.  ?bia  wist  be  done  aiMer- 
Icoliy. 

Vhlle  point  values  of  the  heat  transfer  coefficient  ae  given  in 
Xquaticm  (48)  are  of  acadeaaic  interest,  they  do  not  give  answers  readily 
usable  by  the  designer.  To  predict  the  beat  transfer  characteristics 
of  a  condenser  or  to  enable  rational  design,  it  beccsses  necessary  to 
predict  mean  heat  transfer  coefficients  vhich  can  be  applied  to  the 
condenser  as  a  whole  rather  than  at  Just  one  point,  ivccordingiy  it 
Is  possible  to  define  this  aMn  beat  transfer  coefficient  in  the  same 
manner  as  vms  done  for  the  vholly  laainar  filju. 

(55) 

It  should  be  noted  at  this  point  that  the  substitution  of  q  «  v  h- 
»  r  Y  h-  implicit  in  Equaticsci  (l>5)  neglects  the  effect  of  subcooling 
of  the  condensate  film.  As  has  been  shovn  previously  this  is  very  near- 
ly true  for  the  turbulent  film,  but  it  is  not  tru»  for  the  laminar  film. 
Since,  in  the  present  case,  we  are  considering  both  a  laiainftr  and  a 
turbvLLent  film,  this  will  introd\ice  some  error  in  the  results.  There 
vill  be  no  enror  in  the  laminar  curves  provided  an  equivalent  h*  as 
defined  in  Equation  (10)  is  substituted  for  h-  .  The  error  involved 
in  ncglectin^j  subcooling  may  be  aoeie%^hat  evident  in  the  region  juet 
folloving  traasitloa,  iut  vill  beeane  lass  Soqportant  as  more  and  more 
of  the  fiJja  is  in  turbulent  raotion.  In  Bsost  of  the  cases  encountered 
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tiie  laaiaar  film  and  transition  regicai  will  extend  only  a  few  inches 
from  the  top  of  tlie  plate  and  the  neglecting  of  the  subccKiIing  In  tnis 
Munll  area  ehould  not  appreciably  affect  the  overall  results  fo^  pX&tes 
whoac  lengths  may  approach  several  feet»  It  is  also  felt  that  the  im- 
certainties  of  actxaally  determining  the  transition  point  will  introdiice 
errors  of  at  least  as  igreat  a  magnitude  as  the  subcooling  error,  and 
tterefox*e  this  latter  point  should  not  be  belabored. 

Since  ve  have  aasusoed  a  constant  temperature  difference,  the  de- 
nominator of  Equation  (^3)  asy  be  brokma  into  two  parts* 


u  =  ,.    ,    r  h. 


The  terffi  (t   -  t  )  z^  affiles  to  the  laalnar  portion  of  the  ccaidensate 

o 
film,  and  the  term  (t   -  t  )(L  -  z. )  applies  to  the  turbulent  portion. 

o 
Froa  Equations  (16)  and  (13), 


(57) 
Equation  (^4)  gives  the  expression  for  the  turbulent  length.     Sul><- 
stituting  tb^se  tvo  valties  into  Equation  {%)  yields  the  fixml  equaticxa 
for  the  flUMUB  Jbaat  transfer  coefficient  in  terms  of  functions  of  the 
filitt  thickness. 
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Tbe  Reynolds  {TuDiber  at  the  bottora  of  th«  plate  Is  of  mxut  Interest 
as  It  is  a  dii'ect  measure  of  the  amount  of  v«|>or  that  has  beea  cc^idensed* 


«iis 


K  -pu         and  may  te  readily  evaluated  fx-om  Eq^tions  (^) 


For 
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For  5  <  x^^    <    30 
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n^e       "^  -^        OL  ou  OL.  ^59) 

Ksuatlons  (5^),  (58),  and  (59)  represent  the  results  of  the  tur- 
bulent analysis  in  much  the  scrae  vay  that  £q;uatioas  (l6),  (26),  and  (22) 
represent  the  results  of  the  laeiinar  analysis.    The  same  fonos  of  plot- 
ting are  therefore  used,  the  turbulent  plots  awsrely  being  continuations 
of  the  laminar  curves.    A  description  of  the  n\Baerical  evaluations  of 
these  equations  is  contained  in  the  Appendix  A. 
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B  -  DlSCUBSIOfl  Cy  AMLYTICAL  MBUVm 


The  reauLts  of  the  ttaftly»ls  presented  la  this  tlwsls  are  ccAtaloed 
in  Figures  7  through  12.  'mese  eiirves  are  based  &x  a  rational  approach, 
and  show  the  effect  of  a  vc^r  shear  stress  on  heat  transfer  in  laminar 
and  turbulent  film  ocxidensaticsi  for  a  wide  rani^  of  M^    and  Np^.  Tim 
\i8e  of  these  curves  will  be  lHustrated  later  by  a  aam^le  problem* 

It  is  significant  that  the  plotted  values  of  mean  or  avera^  heat 
transfer  coefficient  include  the  effect  of  a  Inslnar  condensate  film 
on  the  i;is)per  part  c^  the  plate  and  a  turbvilent  fils  on  the  lover  regloDB, 
if  <me  exists.  To  use  these  plots  it  is  not  nacessarj  to  first  inves* 
tlgate  ^idMtther  or  not  turbulent  flov  exists,  this  determination  having 
already  been  suide  in  the  preparation  of  the  plots.  The  values  givea 
are  universally  applicable  and  are  not  d^psiniKit  on  the  geoeetry  of 
the  systea  or  the  fluid  properties  exsept  as  they  miter  into  the  deter - 
Binaticm  of  the  dimensioilees  quantities  needisd  to  «aX«t  the  curves. 

Zn^peotion  of  Figures  7  through  1£  reviaal  several  interesting  points 
concerning  the  aaehaniflaa  of  film  condensatioati .  First,  the  significant 
effect  of  a  vs^r  shear  stare.s  is  apparent.  Tbis  effect  is  more  pro- 
Boiifficed  in  the  lover  IL.  zvnge,  giving  an  increase  in  heat  transfer 
coefficimits  of  several  foM  in  scrvc  eases.  At  hi^har  R^^^  the  effect  is 
leas.  The  effect  on  h  of  the  exact  location  of  the  trKisition  point, 
as  Influenced  by  the  vapos*  shear  stress,  beeaaes  less  si^if leant  as 
the  value  of  IL^  at  the  bottom  of  the  s«urfaee  increaswi  and  is  praetieally 
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insignlf leant  for  le  greater  than  10,000. 

Tbe  effect  of  turbul^ice  in  the  film  is  illustrated  by  the  ehaag* 
in  slope  following  the  transition  points.  It  can  be  seen  that  the  tiu-- 
bulence  effect  on  the  heat  transfer  coefficients  is  greater  for  increas- 
ing IL,  and  for  increasing  Np^.  The  Reynolds  Hueai^er  effect  laight  be 
expected,  as  larger  !L  will  result  in  a  greater  degree  of  turbulence 
and  hence  a  greater  portion  of  neat  transfer  by  the  mechanism  of  tur- 
bulent aixing.  The  Prandtl  Husiber  effect  is  best  seen  by  con^^aring 
the  heat  transfer  plots  for  the  very  low  N   »  0.01,  representative  of 
liquid  metals,  and  the  plot  for  tbe  high  N_  *  10.  The  curves  show  that 
turbulence  has  a  greater  role  in  increasing  the  heat  transfer  character- 
istics for  tbe  hitler  If  .  In  fact  the  low  Prandtl  nwober  curves  for 
the  turbulent  porticai  actually  lie  below  en  extension  of  the   laminar 
curves.  This  may  be  explained  oxice  it  is  realized  that  for  these  liq- 
uids, the  thermal  ccxiduetivity  is  so  large  that  it  becoeies  the  principal 
aeans  of  heat  transfer,  rather  than  turbulent  mixing.  Since  the  aasuBWd 
velocity  distribution  results  in  n  thicker  ccaidensate  fila  than  would 
exist  for  laminar  flow,  it  can  be  seen  that  this  thicker  film  offers 
a  greater  resistance  to  heat  flow  through  cmductioa  than  would  the 
laminar  film. 

It  would  be  desirable  to  cocspare  the  present  an^dysis  directly 
with  experimental  results.  However,  experimental  data  on  filju  cooidetn* 
saticm  in  which  the  vapor  shear  stress  has  been  accurately  determined 
was  not  readily  availsble.  There  does  exist  a  correlatiaa  equation 
baaed  on  considerable  esqperiaental  data,  which  has  alrntdy  bMn  mmtiasMd, 
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•ad  is  giv«a  by  Squation  (^)«  In  order  to  cdopare  tb«  aaauLytlcal 
results  vlth  tills  cc»*reiation  eq^oatic^i,  ttafi  curves  in  Fl^trtm  S,   10, 
aod  X2  hav«  been  cross  plotted  at  a  Heynolds  VUid>er  of  ^000  ftnd  are 
presflxited  as  curves  of  meaa  beat  traosfer  coefficient  versi;^  La,   . 


plots  are  represimted.  by  the  solid  lloes  in  Fi^;^re  1^  for  vario^ui 
Praodtl  Nonbers.  The  correlaticxit  equation  in  the  fona  giveia  by  £^pxa- 
tion  {%)  has  also  been  plotted,  azid  is  shown  by  the  dotted  liiMts.  As 
eaa  be  se«ti,  the  lyeewent  is  qpuite  good  for  N-  «  2  to  ^,  which  lies 
vithin  ti»5  raD|S8  of  2  to  ^  of  the  liquids  actually  tested  and  on  which 
tint   correlation  eqiiation  is  based.  Tbmrt   is  considerable  deviation  out- 
side this  rsa^e^  as  evidesiced  by  the  curves  for  ^^   «  10  and  0*01.  This 
iodieates  that  caution  should  be  employed  iihfta  using  the  correlation 
equation  outside  the  raoge  for  vihioh  it  vas  prc^posed* 

Fl^iare  13  also  gnsphjcally  illustrates  the  eonysiderable  eff eet 
that  a  vapor  shear  stress  may  have  o&  the  beat  transfer  coefficients. 
For  example,  the  h  groiv  for  N-^  «  10  and  'C-  •  50  is  about  1.00,  or 
four  tjjues  the  value  of  0.23  glv«&  for  U  »  0. 

In  the  force  balance  givd  by  Bquatioa  (2)  the  cmly  pressure  grad- 
ient in  the  z  direction  considered  uas  that  due  to  the  height  of  the 
vapor  and  is  represented  by  the  ^  A^  ter&a.  while  this  asau^ptioa 
ws  valid  for  the  case  of  a  flat  plate  in  an  infinite  atmosphere  vhich 
we  considered,  it  may  lose  its  validity  vhen  a  vertical  pipe  is  consid- 
ered. For  pipe  flov,  the  existance  of  a  va^v   shear  stress  is  intimately 
eeonseted  vith  a  pressure  gradient  in  the  2  direction.  To  include  this 
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effect^  it  is  xnem^mmary  to  rewrite  the  force  UiXance^  Equaticaa  (2)  in 
tbe  f Qzm 


r-(i/-iiK^^'^^  ^--"^ 


(60) 


3^ 


To  evaluate  j-^   ,   ccaaaider  a  control  voiuac  containiog  the  vapor 
in  the  core  of  a  vertical  pipe,  as  shovn  in  Figure  1^. 
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FIGIBS  1^ 
Considering  the  vapor  flov  as  one  disiensional  and  inclmiiag  the  change 
in  j^xa^attd,  the  pressure  gradi^it  may  be  evaluated  frooi  a  force  bal- 
ance on  the  control  volume  as. 
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wbere  u     respreients  tht  v^Ksr  velocity. 


(61) 


2 
The  veight  rate  of  flow  of  the  vapor,  "w    «  p    u  it  r .   ,  \rtiich  may 

be  differentiated  to  give, 

The  veight  rate  of  f  lov  of  the  liquid  coiuieiiaate,  v   *  2  «  r  P  vhich 
Buiy  also  be  differentiated  to  give, 

Amc  ^  air  i^o  II 

Since  ^_ii2r  »  -  liiC^  ,  ve  oay  solve  beti#een  theae  t\#o  equatiooa  for  the 
velocity  gradient, 

iii ^  IT  A 

Subatituting  this  ejqpression  in  Equaticxi  (61)  gives  the  pressure  grad* 
lent  as 

(6£) 

Previously,  in  the  case  of  the  flat  plate,  ve  had  taken  the  pres- 
sure gradient  »  (g/g  )  p  •  However,  for  the  case  of  the  vertical  pipe 
condenser  «•  isay  take  an  equivalent  vapor  density 

(65) 
This  equivalent  vapor  density,  -which  talces  into  account  the  presaure 
gradient  for  a  vertical  pipe,  contains  tvc  tenas  wnlch  are  dependGzit 
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on  r.,  the  radiiis  of  the  condsnaate  layer  interface.  For  a  given  tiabe, 
the  equations  presented  in  tne  analysis  can  be  evaliiated  in  a  manner 
•oalogofuyi  to  that  for  the  flat  plate.  SixK^e  -j—     and  u  vill  vary  dovn 
the  length  of  the  tube,  these  must  be  determined  from  the  given  conditicais 
at  each  step  in  the  stepwise  mxaerical  integration*  It  is  apparent 
that  tbie   solutions  so  obtained  vill  apply  to  only  the  particular  geom- 
etry for  vhioh  they  were  calculated  and  for  the  particular  vapor  vel- 
ocity which  vas  ccaisidered.  In  this  procedure,  the  advanta^  of  the 
universality  vhich  vas  obtained  for  the  flat  plate  will  be  lost. 

In  Equation  (63)  the  terms  containing  the  vapor  shear  stress  and 
the  vapor  velocity  are  subtractive  and  my  quite  po&slbly  largely  can- 
cel out.  If  the  net  magnitude  of  Equaticn  (63)  is  snail  eoDpared  to 
the  condensate  dmxsity,  p,  it  is  then  possible  to  apply  the  results 
obtained  for  the  flat  plate. 

It  Is  ^forthvhile  to  investigate  the  magnitude  of  the  pressure 
gradient  correct icai  for  a  specific  case.  Consider  a  1  inch  diataeter 
tube,  in  vhich  lAter  is  coEUiensing.  S&tiirated  vapor  teaperature  vill 
be  taken  %b   120  F.,  and  the  tube  vail  t«Bqperature  as  60^  F.  To  take 
an  extreoe  case,  assusoe  "Hr  «■  ^o. 

Itean  film  to^erature,  froa  Equation  (1), 

t^  «  120  -  3/^  (60)  »  75*^  F. 
Then 

P  -  1.000  X  10'^  ft^sec 

vkj  *  ^'^^  X  10  ft. 

Since  /^  <'</ 
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Frooa  Equation  (^2), 

UL  ,  J_  /t)  =  _L  ii  At 

ii) 


w- 


Consider  a  point  in  the  tube  where  x^  ^  ^0»     Then  x   b  500  and 
0,1 


At  75°  f»,   the  Prandtl  !Jiaaber  of  water  ia  6.29,  Proceeding  in  the  aan* 
Gter  outlined  in  the  secticoi  on  calculation  methods,  F^  laAy  hs  evaluated < 

Fg  •  57.6 
!Ihfin  taxing  k  of  the  condensate  as  ,557  BTU/hr-ft  -T, 

h  «  2660  BTU/lir  -ft^  -  °F. 
Substituting  h,  h^  «  1026  BTU/lb,  and  ^t  «  (yO^T., 

-V-  «  2.59  Ib/hr  -  ft^ 
Taking  «BOOth  pipe  friction  factors  for  illustrative  purposes , 
f  »  .005  at  a  vapor  Reynold's  Number  of  64,500.  Since  Tat  =  t  ^"t  - 


v»  say  solve  for  the  vapor  velocity,  u, 
u  ">  1090  ft/sec 

^^e^v*  ^^^      =»  64,500 
Substituting  our  oon^mted  values  in  Equation  (65 }»  we  are  aible  to 

evaluate  p  , 

P„*  »  oi"»  >  2.^  -  21.32  »  -  19.4a  lb/ft ^ 
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Although  p  •  i«  about  l/5  odf  th©  condensate  density,  p,  thia  la 

not  tfa0  BMg»ituA»  of  tbe  error  involved  hy  neglecting  this  term.  Be- 

feiTing  to  Equation  (oO),  by  continuing  our  i^reaeot  specific  example, 

we  can  evaltmte  ti»  maximua  error  in  the  sbear  stre»«,  T  ,  which  will 

oceiir  when  x  «  0,  or  T  =  U  ,  since  x   »  50,  x  «  50  X  1  Jwj  x  10 
'  o     '  o 

ft.  «  .0073  ft. 


To  »  .^55  +  (62.^  +  19.^)  (.0073) 

Cyc,^.'  .^55  ~    .597  «  1.052  Ib/ft^ 
By  neglecting  tbe  pressiure  gradient^ 

(ToJ  ,.co.n  •  -^55  *   .^55  -  .910  Ib/ft^ 

Sbe  final  error  idliich  is  introduced  in  this  specific  case  is  of 
the  order  of  15.5%.  These  calculatiraas  are  inclined,  not  to  exhaust- 
ively Investigate  the  e^ect  of  t±»  pressure  gradient  term,  but  rather 
to  give  an  order  of  asa^itude  of  the  error  vhich  might  be  expected. 
The  values  chosen  in  the  probl«EX  ai:*e  felt  to  be  i-epresentative  except 
for  the  vapor  velocity  of  1090  ft/sec.  Iliis  high  vaiiie  vas  deliberately 
picked  to  indicate  aa  v^ppcsr  limit  to  the  error  involved*  In  ajoy  prac- 
tical steam  oondmoser,  the  vapor  velocity  will  undoubtedly  be  lauch  lower, 
with  lower  vapor  shear  stresses.  It  may  be  concluded  that  in  many  cases, 
the  eiYect  of  the  preasure  gradient  aiay  be  neglected,  and  the   analyti- 
cal results  for  the  flat  plate  siaj-  be  applied  to  vertical  tube  conden- 
sers with  cenfidence. 

m  all  of  the  calculations  for  Figures  7  throu^  12,  the  quantity 
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(X  **  Py/p)  ^<^9«  tak«i  M»  uaityj  <8«g*^  p  /p  floall*     It  vas  also 

ttetUand(t      .  t  )  <«»  u»ltom  alcne  tt»  plat«.    -m.  •v^tioo. 

o 
la  the  aaml^ls  mjb^  b«  r«Mi&lly  €¥mImt»A  aoecmatiJEkg  for  any  value  c^ 

(l  -  p^p)  and  for  aay  p£r«»crlto<ftd  varlatioea  of  (t       •  t  )  or  of    Ta^ 

o 
alockg  the  I<mgth.    Tbo  affacts  of  tbeae  varlaticMa  ara  fanctiocia  of 

tiu6  particular  etetmtry  and  ara  hard  to  ganaralisa. 

Ot  the  mfiiay  aaaueq^tlooa  voderlyisg  tMa  anal^rals,  aoat  ara  mamiA 
azid  have  cooelderahla  axperi»«ntal  avidaoee  to  back  ^£>  their  validity. 
The  tvo  most  subjeet  to  quest  ic»ri  are  tSam  aaaveed  valocity  dlatributioxx 
aod  the  criteria  for  traaifiltloa.     It  will  be  t^te  purpose  of  tim  eecctid 
aactlon  of  this  theals  to  l»vaatlgate  the  bydrodyaaailcs  of  a  falling 
vertical  liquid  fila  iBider  the  influftace  of  an  a^^praci labia  ahaar  stxfsss 
at  the  interface.    Ttnnm  iBvaati0atioo«  ara  lurallaiaary  oaly  aad  vlll 
Bcraly  point  out  the  directioa  to  be  follovad  in  a  aora  detailed  exper- 
imental progjtn. 


-0^- 
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Application  of  tisM>  r«raltft  of  the  an&lysls  is  eonvenleatly 

vith  the  \w«  of  tlis  disa^naiosiless  plots  of  Figures  7  through 
X2.     Sixuse  aean  or  ovorall  IooaX  traa«fer  coefficienta  are  daslt 
vlth  throi^hout,  ci^  point  by  point  or  9tepvi8«  ealeulAtioas  9r« 
necftssary.    Tba  mstsml  methoda  involvHed  in  a  design  calculation 
are  illustrated  by  the  following  probleas* 

The  tvo  problems  seleeted  are  felt  to  be  representative  of 
those  vbich  slight  be  ast  in  industrial  prectice .     The  operating 
eonditioBA  and  tobe  slsee  chosen  eorrespond  to  two  actual  eases 
mMnrs  SMNhMBTSnsnts  have  been  talccna.     The  caleulaticms  should 
therefore  serve  the  additional  pitrpose  of  eooiMurlng  the  aceuraey 
of  the  predictions  asds  by  the  presMit  analysis.     The  first  prob- 
lem utilizes  the  initial  data  for  \eater  run  I9I  of  Carpenter's 
•spwrlwRital  progrsBi,  reference  (3)«    The  seeoc»i  is  based  on  an 
actual  Indus  trial  desi^  as  reported  by  V»\Krster,  reference  (1^). 

The  problen  is  posed  in  the  following  amaasr*    Tor  a  givm 
condensate  iieight  rate  of  flow  per  tuJje,  it  is  desired  to  deter- 
mine the  required  length  of  tube  to  ccnidbsasft  the  amount  speci- 
fied, vith  no  uQeca3d<msed  vttjpor  passing  out  of  the  bottom  of  the 
tdbe.     In  addition,  it  is  required  to  compute  tba  »••&  enrvftwising 
side  h(»t  transfer  ooefflciects. 

Jn  both  cases  the  ccaidensing  vapor  is  steam  at  essentially 
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atmospherlc  pr«ssiire*  The  wall  temperature  has  been  takma  as 
the  arithmetic  average  of  the  wall  teoiperatures  reported  In 
each  of  the  references. 

So.  1       No.  2 


^bc  Character i«tle0:  >/8",l4  BWG   l',17BW 

laaide  Diameter  .i^59        *^^  in* 

,  -       2 
Craw -Sectional  Area  .00115      .00426     ft 

Tube  wall  tcnperature,  t^  162°F     182.8°  P 

Desired  condensate  flow  rate,  v         >0.1       201       Ib/hr 

—  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  _  —  —  —  —    p 

Vapor  laass  velocity  at  top  of  tube,G.  ^5,600      ^7,100     Ib/hr-ft 

Since  Gg'O,  Eq.  (57)  gives  G^^       2?, 300      27,^00     Ib/hr-ft^ 

{tpseific  voluroe  of  saturated  vapor  at  l4.7  psia,  v«26.8  ft  /lb 

Mean  vapor  velocity,  U^,  186.2       20^      ft/sec 

Cendeosate  flow  rate  per  unit  periiaeter 

at  bottom  of  tube,  f  -wA  ^l3       870      Ib/hr-ft 

Average  condensate  flow  rate  Vt<^  209        ^55      lb/hr«ft 

Since  vater  is  the  fluid  being  condeaaad,  the  ratio  of  surface 
tensions,  s  /s  =  1. 

The  density  of  saturated  liquid  at  212°  P.,  p  =  59*7  lb/ft' 

Then  ^  5-5        7-^     ft^/hr 

Tl»  viscosity  of  st««a  at  212°  P.,  ^  »  .0303  lb/ft Htir 

The  Beynolds  Nuaber  of  the  vapor,  -j^       29,700     66,600 

Utilizing  the  tvo  ccmponent  flov  friction  factor  data  in  reference 
(U),  f  ^     .0063       .0072 

Vapor  shear  stress  C^~   t       g~~     .1296       .17^     lb/ft 
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mva.  Film  Temperature,  from  Eq.  (1),  t^      ijk,'}^  190° 

Kinaaatic  vi»co6ity  at  t„,  /^  .599  X  lo"^  .562  X  10*^  ft^ec 


f 


jH:)^^^  •T90  X  lo"*'*  .7^1  X  lo"^  ft 


.^_  Sot. 


Vlacoaity  of  ccaidencate  at  t^,  JJ  .75^  X  lO'^'   .680  X  lo"-*     lb-«©c/ft' 

Condensate  Reynold*  Jhitjer  at 

bottoB  of  tube,  Np^  -  ^  1915  ^20 

Prac^tl  ^t»ber  of  Ccsidenaate  at  t-,  K-^  S.26  2.C© 

Enter  Fig.  10  for  li^   «  1  vith  H^^ 

•nd  r:   ,  get  Km  ^-£)'^^  .^1  .^7 

Enter  Fig.  12  for  ¥.^   »  10  with  H^^ 

and  rv  ,  get  if  {^f"  .7^  .98 

Interpolate  to  actual  Sj^,   b^ /'-^  ^  .^  .55 

Itoermal  Conductivity  of  condenaate  g^ 

at  t^,  k  .386  .390   u,,.ff°F 

Ntan  heat  transfer  coefficient,  h^^  2250  2^70 

solve  for  z            ^  5.60  10.4     ft 

Actual  z,   as  determined  by  experlaent, 

(references  (5)  and  (l^)  3.6  11.7     ft 


By  cc«q?iarios  calculated  and  actual  tube  lengths,  it  can  be 
that  "Uae  analytical  results  will  p2-edict  the  perfonaance  of  a  coiKleaser 
quite  cloaely.  ProbaUly  the  greatest  source  of  error  lies  in  the  selec- 
tion of  the  coarrcsct  average  vapor  sheer 
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8tr«ss  ei3d  tube  vkU  toi^rctture.     If  it  vas  d«»ired  to  predict 
perfona&nce  vithin  very  close  tolerances,  tbe  mtsmrieml  inte* 
gratlon  preiMMtant  described  in  Ai^ndix  A  could  be  reworked  vith 
the  Imowa  vsriatiosis  of  ^  t  said   ^'^   inclufded  in  a  stepwise  a^B^gnnM" 
Smaticax,    Bovever,  this  vouM  involve  a  rather  tediom  end  lengthy 
set  cf  calculations  and,  for  siost  eases,  it  is  felt  that  tJas:  cosi- 
stant  average  valiies  as  iised  above  vill  porovide  resiii^  of  suf- 
ficient accuracy. 
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PART  TI 
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PAHt  II 

Beurt  I  of  this  theslii  I0  ccuacemed  vith  tiie  development  of  aaa-> 
lytic&I  expresaicKUi  predietiag  tjae  beat  treosfer  coefflclente  «xf  a 
fallijD^  film  coateBMr.  Use  has  beaa  oade  of  the  heat  lacnea^ttn  trane* 
fer  aoology^,  so  that  the  validity  of  the  predicted  heat  traasfer  coef « 
f  ici^atfi  is  dependant  upon  a  looowledge  of  the&  h^r^'odj/nauice  c^  a  falling 
fills.  AMmmfftimm  have  beea  aade  regarding  the  behavior  of  such  a  f ilA, 
based  on  the  beat  Jmovladge  available  toda^.  It  is  the  purpose  of 
this  section,  ti^refore,  to  give  a  brief  reviev  of  the  literature  cov- 
ering the  subject  of  falling  liqiiid  filias. 

Moerad  and  Badger  (1^)  diseuss  the  onset  of  turbulence  in  a  falling 
liquid  fils.  Their  cCM:u:Ii4Sion,  based  on  the  work  of  other  investigators, 
is  that  the  critical  Reynolds  Hlflber  is  in  the  rso^  of  l300.  They 
also  cite  data  of  Jakob  and  Erie  (19)  sa  shoving  that  transition  osjr 
occur  at  lover  21.  in  the  j^resenee  of  high  values  of  vq^t  velocity. 

In  reference  (IT)^  Eebbard  and  Bad^gsr  aaaaared  film  coeff ieieaats 
of  heat  transfer  on  a  1  inch  O.i).  vertical  tube  condei^ser  12  feet  long. 
Beynolds  B««^»ers  of  the  filiE  covered  the  ran©e  500  to  !iCX)0.  Results 
are  consistmitly  higher  than  thoiw  predicted  by  the  Nusselt  Irimlaar 
efuations,  and  gwoerally  lie  in  a  line  parallel  to  the  predicted  lijoe. 
Ho  Bsteaaureeaants  of  the  rupor  shear  stress  vere  mmSm,  md  it  is  felt 
by  this  euthor  that  the  discrepancies  may  quite  possibly  be  eKplained 
by  this  effect.  No  infonsiation  as  to  the  ctsset  cf  turbulent  flow  osa 
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b«  dcdi»5«d  frcm  thes«  results. 

KirUbrlde^  reference  (2),  jwreienta  heat  tracefer  data  of  several 
Inwstlgfttoani  cmq  liquid  filsuB  flovlag  do^m  vertical  \ealls.  This  data 
eovere  the  Beyoolda   msaber  range  of  700  to  20,000.  By  tinB  narlced  dif- 
ference In  slope  of  the  correlation  line  fron  the  HUssalt  liwinar  pre- 
diction, turbulent  flov  of  the  filtt  is  definitely  indicated*  Hc\«ver, 
there  is  little  data  for  the  lover  I^^,  so  it  is  not  possible  to  airect- 
ly  Infer  the  location  of  the  transition  point*  The  best  correlating 
line  through  the  data  intersects  the  theoretical  Ihi^aelt  predieticm 
at  a  U^    of  eboat  1100.  Tim  vapor  shear  stress  ism  not  coosidsred* 

In  addition  to  considering  beat  transfer  coefficients,  Kirkbrids 
(2)  presents  data  on  Bsasured  fiLa  thictoeseei  tw  IL  less  than  £000. 
Excellent  agyawftwat  vith  predicted  thicknesses  is  obtalzked  for  the  lover 
Nj.  ,  hoveirer,  deviations  are  observed  for  the  higher  8^,  .  Surface  iMves 
are  given  as  the  esqplanation  for  the  dlserepancies. 

Cooper,  E^rev,  «iif&  MBJkdsns  (iB)  present  isothermal  data  o£  six 
observers  dealing  with  liquids  flowing  down  the  wetted  ssaooth  sun^nee 
of  vertical  towers  or  over  flat  plates  at  sli^t  inclinations.  Fil« 
thicknesses  vere  seasured  either  by  a  aderometer,  visually  observing 
contact  vith  the  f iljs,  or  by  suddenly  stopping  the  flov  and  aeasuriag 
the  fiisntity  that  drains  out  of  the  test  section.  In  a  plot  of  the 
Fanning  friction  factor  versus  n^  .   scatter  of  th«  data  from  the  lami- 
nar equation  starts  at  a  5L  of  about  1000,  with  the  best  correlating 
curve  indicating  transition  at  about  2100.  Data  ^«re  not  available 
for  determining  the  i^Misible  effects  of  ^fOM  velocity  aa  the  liquid  flov. 
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Conaideretole  exptrisaentaX  vorK  has  heen  done  on  non-condenaing 
lieat  exchangers  ttnploying  vertical  tubes ,  in  vhich  the  fluid  flovs  in 
filn  form  dowu  the  tube  valls,  refereaces  (20)  to  (26).  Bays  in  ref- 
erence {2'£)  has  plotted  heat  transfer  coefficients  versus  N^  •  A  break 
in  the  curve  at  a  N-  of  1300  or  l^MX)  laay  possibly  indicate  trans  it  icm 
of  the  film  from  laminnr  to  turbulent  flow,  bovever.  Bays  :i  sriiWHiiids 
a  value  of  1000.  Leavitt's  data  (25)  for  vater  extends  froB  a  li   of 
1000  up,  with  no  Indication  of  transition  as  lov  as  1000,  although 

eoosiderable  scatter  of  the  data  is  observe  belov  N^  •  1700.  The 

t^e 

data  of  Fentress  {2k)   covers  the  ReyiK>lds  WmBbw  range  of  350  to  20,000, 
vith  no  break  observed  which  might  indicate  trans it ic»i.  In  all  of  these 
references,  no  vapor  shear  stress  vas  present  and  no  atteispts  vere  Biade 
to  aBsasure  filen  thicknesses. 

Dukler  in  reference  (27)  has  used  a  film  thickness  meastaring  tech- 
nique eiaploying  a  capacitance  neasureaient .  This  device  determines  film 
thickness  by  the  diffeirence  in  dielectric  constants  of  air  and  water 
beti#e^i  tvo  plates  of  a  capacitor.  One  of  the  advantages  of  the  tech- 
nique is  that  it  averages  out  surface  wares  and  irregularities.  All 
neasursments  viere  madis  on  the  filta  flowing  dovn  a  vertical,  flat  plate. 
Film  thickness  readiasa  affree  vail  vith  the  Huaaelt  predicticm  in  the 
lower  N   range  and  vith  the  prediction  from  the  "universal  velocity 
distribution"  in  the  higher  H   ranse.  As  t!ie  tcKshniq^e  la  aoobawtiat 
difficult  to  6$ply  in  a  round  tube,  no  investigations  of  a  vapor  shear 
stress  effect  vex*e  laadertaken.  For  the  flat  plate  with  no  veq^or  shear 
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strees,  traaaitiaa  Is  indievtwd  at  a  !i_  of  1000  • 

It  is  a^^parezit  timt  ccxuilderablft  variation  «xi«t«  la  tbe  results 
of  the  nany  i&vestli^itora  elt«d  b«««    On  tba  basis  of  tlie  v»rks  cov> 
•red  la  this  survey,  oo  dlf  laite  eonelusloas  eaa  tm  nwehed  as  to  tiie 
oameit  of  turt>uUttee  in  a  liquid  f iljoi  f loving  dowavsrd  ui^Uur  tbe  influ- 
ence of  gravity.    FurtbenBore.^  no  Invest i|gatl(»s  faanre  bean  bmuIs  to  de« 
teralzMs  tlie  effects  of  an  Bmprecisble  ra^iop  i^bsar  streas*    Hoare  kaowl* 
•dga  in  this  direct ioa  is  urgeatly  needed*    This,  therefore,  vill  he 
the  direetioa  to  be  talwa  i&  the  eaqperloMiBtal  part  of  tha  thesis. 
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The  puzpoMBi  of  the  esq^ttrSitumtflLl  program  is  to  invvstji^te  the 
h^drodyroBLiea  ot  »  liquid  fllsi  flxrvlng  donnward  wHooq  a  vertical  our* 
fae«  VBOder  the  Inf luisiee  of  a  sbttar  str«M  extrted  on  tJM  cmter  f ao« 
of  l^e  liquid.     In  order  to  obtain  queintltative  reeults  it  1«  imnmmb* 
sary  tiQ  iw  able  to  directly  measur'e  the  foIXowii^  fuantitias  $ 

(1)  The  liquid  flow  nate 

(2)  The  li<2Uid  tniperaturtt 

(^)    'Th*  awn^itude  of  the  lnt«rf«eial  skiMkr  str«s« 

(4)  The  Uquld  fil»  thie)UM»0 

(5)  In  addition,  it  voolA  b«  dissixii^e     to  b«  able  to  obtain 
sMMumreBWjata  ^ieln  >«uld  |i;iv«  fUMntitatlva  infonafietloa,  aa 
to  friction  factors  for  t^»o  ceagaaattt  flov 

Although  a  tvo  dlBMBaBlOBal  aet-up,  vith  t^e  liquid  flowiag  doim 
a  flat  plate,  mm  eonaiaerad,  the  aaeaaaity  of  accurately  determloiog 
the  interfaclal  shear  atreas  lad  to  the  choice  of  a  vertical  tub«  for 
tha  tast  aaotion.    By  introducing  tlia  liquid  along  the  tiA>e  wall,  it 
flOKa  dovB  the  valla  in  film  form.    A  flow  of  aiz*  in  the  coo?®  of  the 
turiM  aiitKULates  tha  m^por  f  lov  in  a  vartioal  tube  condcnaar  and  sillies 
the  ahear  straas  on  the  outer  face  of  the  liquid.    By  aaavniag  ao 
preaeure  drop  across  the  liquid  f  iln,  t^iia  ahaar  atrasa  cart  ha  eaaily 
WHWured  by  pr^iam'c  drop  r«adin£p»  along  the  Itmgth  of  the  tube. 

1^  t^t  section  is  shoiai  in  Figure  13  an>^  la  a  c;ne  Ixich  diaoeter 
tuba  approaclBMiiteXy  '^  feet  long*     Pressure  tapa  ar@  installed  at  one 
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a«>xiq<suc  Sine  '\mmjaaiiSU3io  e^oua  laiUiTimi  «  iit  voi^  -o^sv  0^7  tm^asUamin  9au9 
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foot  intervaXs  olons  it«  length.     In  oi-der  to  prevent  arcing  over  of 
the  liqixld  at  hi^  tlont  rates,  a  tnooth  coavm^pmt  noszle  is  placed 
at  tlie  top  of  thtt  tube.     Surrounding  t^ie  naizzLa  Is  a  settling  tank 
tc  castsp  out  ax^  turtouleruse    or  f  luctiaatlatis  ¥hich  sili^t  exist  In  the 
liquid  as  it  is  dlselaarged  froa  tHe  sugpiidy  systos.    !tbe  nozsle  acts 
tm  a  veir  In  the  settling  tank  to  flMilS^  distribute  tl^  liquid  along 
tte  vaXls  of  xaa  tul>e»    TiM  entire  test  sectiou,  aozzi»,  and  settllag 
tank  are  eoostructed  of  lucite  to  permit  visual  ^wervation  of  fli>w 
oliflQraeteristles  • 

Six  probe  ports  are  iastalled  in  the  test  siKstlon.    These  are 
arrsaeed  in  two  groiq^  of  three  each,  spaced  evsaily  around  the  cir- 
ei»fer«nee  of  the  tuiM».     Details  of  the  probe  ports  are  shoiiai  in 
Fiipure  16.    Flla  thl&kneas  readings  are  made  with  the  oodlfied  depth 
tms^  miaratmtw,  shovn  in  place.    Metal  probe  plugs,  "^^joae  inner  ends 
are  flush  vith  the  test  seetion  ^mll,  are  installed  to  provide  elec- 
trical contir^uity  from  the  mlffmastsr  probe  through  tOae  liquid  f  11a* 

The  aieroB»ter  probe  is  a  standard  depth  ga^se  atlerosieter  vith  a 
pointed  iractensicxi  added  xo  'tihe  shaft.     The  shaft  and  the  exrteosion, 
except  for  the  cxtresae  tip,  are  insulated.     This  probe  is  Bhowx  in 
Figure  17*    To  obtain  filA  thieioieee  readlxt^^  the  probe  is  Inserted 
in  tlk  port,  vith  the  sboaliiers  of  the  iiii^v<Qaftter  held  f  iiialy  on  the 
blocks  proridsd.    A  sero  resii.ii^  is  eibtained  by  advancing  the  micro- 
■eter  until  it  touches  the  fhr  vtall  of  the  test  seetlcxa.    The  pxx^be 
la  retracted  boA.  a  iW(isitiV«y  jOlBMter  cimnected  between  the  end  of 
the  Qicrotaeter  aad  one  of  the  adjacent  port  plugs,    the  mici-oaeter 
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tettilq  «i  ttiaaoa  ^jngnrvoM  4Uoai»  «  (••to'i  v  J- 

iosi-  aoli^iMi  «  ml  •laimi  Mtt  aaiJksiovxve    .  9tli  tm 

ma  al  t»lM»  SOi^lm  ttakdw  m»l$MuS'*'-*''^  ^  o* 

eJ-M  tXjMOM  «dr    .flKNPia  VtiW  *^  ^^^''^  '  aJt  11  m  sx 

aoouU  lilivli  Mtl  Aiiitftv^v^t  \lmft»  cS  ^s^  al  tlmi 

•  9ol^%l^M^mnmi» 
919  •mmO    .aoilMa  i9mS  mu  al  b9U*^%al  mm  a^soq  Mfvaqi  xtl 
-119  fldi  Smucrm  xLmtfn  toMqpi  ,dA»9  (!>viiit  1x>  •VKyii  curl  Al 

tfi^  teltitos  mU  tfllv  •&•«  wrm  ■nfflAMTx  (MMotaitSl  «£1^    .51  •va^l 

«te»  iSAOl  Me^  ««y(Xq  9d&tq  l»tm    .•&&£%  oi  iaf€^  ^'mttmyt^tm  t(ym 

^ftoiBmsh»  iMfl  bna  ItaSi  «it7    .nMlft  «Eftt  ol  teAito  oslMMJm  b»$aloq 

teluttil  al  94^n  mdS  ^obHww  ■■WMfiirrtt  «Ut  ftkaltfo  oV    .Ti  ni«ll 
•cil  no  xlsr'it  bi'»€  t»&imn^m  Mil  V>  r»&£irail»  Mil  dilJtir  «lwq  wdt  al 

-O'l-u.    ^■^0    ^iLT'Ti  .;v.,t;    X*  6««l*l«fe  •!    »«**••*  OWHI  A      .bmttV&Vl  93LS>Old 

»(fc-:q  i»rfr    .ftoi:,ti>»«  HmT  •rf*  'to  XXfiir  iftl  Mil  Mdanol  11  Illaw  imlMi 
"i<;  bna  .^rt,:^  r*?njtj?^  b»SrmmMH>  ^vliMwuft  tTTtiiaiwi  «  Ims  telfMrchn  al 
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iB  then  advatMsed.  tiatil  eleeti^ioal  continuity  is  emtebXiBtysd  by  the 
prc^«  tip  toueJiiag  tJMft  outer  ed^  of  tbe  liquid  film.     Thia  contin- 
uity i»  aasily  dotectcd  to^  the  ohometar.     The  Insulation  cm  the  shaft 
of  ths  yrolM  pr^virntft  cocapl«tioa  of  the  circuit  l^y  the  f  Ua  oc  the 
BMur  uall  of  the  tuhe^  through  vhieh  the  parabo  extends* 

pyeisme  reediags  in  the  test  secticu  ^M»r«  made  with  a  stsndsgd 
U  tube  mancoieter  with  air  in  the  tubing  betwcea  the  stfuaopetor  «nd  the 
presMre  tap.     Surface  tension  forces  at  the  interface  betweesa  the 
Xlquid  tm.  the  vails  of  the  test  section  and  the  air  In  the  tubing, 
plus  tlw  twsAancy  of  the  liquid  to  be  dravn  froR  the  test  eectlon  in- 
to the  tuhln^  utun  sJMSiglng  ttee  ppensure  Xevel,  oade  tiiei*e  readirD^i^ 
■MBB'wfaatt  difficult.     UDcertaJiities  of  waasured  pressvire  levela  vex« 
oalj  evident  at  lov  air  flow  rates ,  lAiea  the  surface  tension  forces 

r^ifsre  of  a  ccisparabLs  magnitude  to  the  pressure  forces.  Modiflo«44ani 
tere  been  nads  to  correct  thla  by  using  a  siaglje  mmmmimr  tube,  vith 
the  eonoeetiBg  'teiblag  coupjuetely  filled  with  Xiixald.     This  requireo 

vMMKt  s  sero  reading  be  taken  ulth  no  flcv  in  th«  test  secticfi,  but 
ellsilnstes  all  interfaces  in  tine  pressure  neasurlog  opparatubB.     Hfow- 
erer,  all  data  9smmt^mA.  hare  uare  aade  vith  th^  orl^i^lael  nywteR. 
Pressure  readings  plotted  as  stanklght  lines  for  all  but  the  rtssy  low 
air  flow  rates,  and  snpeai*  to  be  fairly  reliable. 

Hkm  test  secticm  diaeharges  to  a  olosad  aaparatlng  tank,  whc^e 
purpose  is  to  separate  the  liquid  from  the  air.     The  air  suction  line 
is  connected  to  the  top  of  the  tank  through  a  standard  pipe  Tee  Inside 
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hiiafcwf  a  ditv  •bmm  avsv  aoitoM  ^m^  «ii  ai  npnll— i  rsMMtf 

Mi;r  xUMw^acf  •ouClis^aJt  9dS  $m  •Mrvcd  nolaoKT  mm^iu^     .qpitf  niNMiaf 

«yilclki#  «1#  xUt  xta  1^  Am  ooi^Ma  Smit  tdst  lo  ail«v  mI^  jk>  Aiayitl 

-at  attA;r»<ka  :»«i»^  Mf#  a^^lt  tfnn^  ad  o^  himptl  wtU  \o  x>ftii«nt  «d^  mtlq 

i»nfrh— 1  oAMl^  ttftMi  ^KTal  •waaanDi  ail#  fli— ill  aaAr  yitda^  atf#  ot 

•w  alonral  wmaan^r  bmmmmMi  Jo  aftitixtjKhtMMtf    .^iiir9llli6  ^artwaaai 

aaovot  aolaoai  aaaftwi  adtf  eacfir  ^ad^Tjen  noil  xta  vol  ^«  Simttn  xLtn 

mmk*M9XMtxM    .aaamt  avnarsg  adl  cuT  ate^Zaq^Mi  •Iilata^pin  a  to  Mtir 

ilc^Jtv  ,MlM^  falaaBiifli  %I^^%  m  j^lau  ^  aiil^  Stnmcf  o&  •Ami  aaMf  arad 

•m.tup9r  •IdT     .btapU.  d#tv  MXXlt  x4a^«£^pr>9  a^dk(#  Boltoaoaoo  a4^ 

j'Mf  ^aoUvmm  laai^  aatf  «i  volt  oo  M0kw  ouiMS^  atf  solteait  <n»a  a  AaAT 

-voK     .avte-vavva  iat'suaaa*  araiaaaai  ««ft  cU  aAa«t*ui»ai:  ixm  —imaimtlm 

,mtwff  Xaolal-Go  ad^  dtkn  ateN  ••sair  a^»w.  .^^ir^aattattq  mtMh  Um  ^lara 

val  tiaT  ad#  ia^  llm  -sot  aaeil  m^JtmiS*  em  bm&^olq  u^tkmn  aiaaa'rt 

•  •ic  '  9d  ot  ymaf^pi  tea  .aadarx  troll  oJU 

aaoi^v  ^Mamif  ;$al^»'mtt»m  l>ft«v  ^^^  4,  >^t-  aa«M»*^j.j.>  lAoid-oaa  ta*#  adf 
aolX  acid-au«  iU  aeff    .xta  adt  asc^tl  :s>tupil  94i  aAaiuHpa  <^  ai  aao^paiq 
a£itiai&t  aafl?  a^.  J^awi-a  a  d^iponni^  jdat^  adt  Jo  qo^  mtf  e#  .hairaaw^a  %l 
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the  tank.     The  air  flov  oiust  traverse  several  sx^arp  chat^i^s  ir.  direction, 

thereby  separatii^ifi  miy  liquid  droplets  vi^ieb  mln^t  be  entrained.    A 

liquid  recirculating  pus^  is  co;meoted  to  tiM  bottom  of  ti:ie  separating 

tauk  ti-jToiagn  a  flexible  iioee. 

Mo  atteiapt  was  made  to  measnxre  t.ie  air  flov  rates  ijn  Xbmse  prelim- 
inary investigatipMS .     Sue  a  aeasur^ieuts  could  be  easily  made  by  instal'* 
ling  a  rotaiBeter  or  some  positive  displacecieat  ateasurii^  device  la  the 
air  siaetioa  line*     It  is  felt  tibat  tbese  oeasureaents,  Wiich  vould  ea* 
able  deteraiaatioti  of  friction  factors,  vould  be  aigaly  desirable    for 
future  experinaeutal  work. 

Tkm  liquid  flow  rate  is  measured  by  an  orifice  installed  in  txie 

siqpply  line.     Tint  pressxare  drop  across  this  orifice  is  Bteasured  by  a 

i 
water  over  mereur;^  laanoowter  •     This  systsa  was  calibrated  in  place  by ' 

tiae  "  wsi^it  afsa^ureoents .     Calibration  data  aiid  curves  are  given  in 

Aii^ndix  B.     Liquid  tei^peratxure  readicgs  are  obtained  froa  a  tbemocaeter 

placed  in  tJae  settlioi^  tarJc  at  the  top  of  the  test  section. 

Auxiliary  equipnent  consists  of  a  pa^  to  recirculate  the  liquid 

through  the  systsA,  and  rtecessary  piping.    For  the  preliniriary  experi- 

lae&tal  progrsn  initladsd  in  this  thesis,  water  was  the  only  liquid  tested 

arid  therefore,  connections  \mtt  made  to  a  vater  laain  as  a  soarce.    A 

line  sketch  of  the  equipoent  is  shovn  in  Figure  l8.    Figure  19  is  a 

photograph  of  the  set-uqp,  nbaiinjs^ iSHt^mt  seetioa  assesibly  ar«d  sepa- 

ratija^  tank  oti  the  ri|sbt,  with  the  micraaeter  probe  in  place. 


J 
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Since  it  was  aatieip&ted  that  vibrations  vcmid  quite  possibly  be 
critical  in  affecting  the  steadiness  of  the  flov,  flexible  hoaee  ifex>e 
installed  in  the  pipi^^^  from  the  separating  tank  to  the  pun^  and  in 
yie  Buqpply  piping  Just  before  the   settling  tank.  TtM   test  section  was 
thus  isolated  from  vibrations  In  the  piping  system  acid  puap* 
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IHETION  C   -  FI04  THICICSESS 


riXni  thickniBftS  wBsm&rmmniM  have  been  naOxs  in  the  maxoBT  dyeecrlbed 
lA  SMStioB  By  Part  II.     .Since  ve  are  dealing  vitb  a  x-cnmd  tube  rather 
than  a  flat  plate,  it  vould  r>e  worthwhile  to  inv««ti^te  the  •q.iMbtlOQS 
iqppllcabiG  to  this  case  and  coaqpare  tV^eaai  vith  tiaos^  derived  for  the 
flat  plate  In  earlier  sections. 

Consider  the  control  voline,  ithoae  cross  action  is  aho^m  ahadad 
in  Figure  20*     This  control  volume  is  an  aanular  ring  aroiaid  the  tube. 
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FIGURE  20 
Since  there  is  no  coodensatlon,  maaentua  tenaa  ma^  be  neglected, 
and  the  force  balance  for  the  control  voluae  is  given  as 
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(64) 
After  ftocBO  •Iffipllfication  and  rearraiig«sa«nt,  this  reduces  to 

(65) 

Since  moiaentua  tenas  ao  noi  enter,  it  is  poeaible  to  use  the  eoc- 
presBion  for  ^^/U  of  Equation  (61)  vith  the  secorid  tcnc  «  0 


(66) 


I%ft>dtitutii3g  this  in  E9as;tiQa  (65),  v»  gst 

(67) 

3ii  solving  for  dp/dz  from  cooditicaas  in  the  vaijor  core  of  the  txSbe 
SBd  thMi  applying  this  dp/dz  to  the  equations  for  slimr  stress  in  the 
liquid  film,  the  aasiavtion  has  been  aade  that  ao  change  in  the  pres- 
sure sradlent  occurs  in  a  radial  directioa.  This  assisa^-tioQ  is  ^ite 
good,  and  alBiOst  exact  fcr  thin  liquid  filios. 

In  order  to  compare  the  distribution  of  shear  strees  given  by 
I^uation  (67)  vritli  that  previously  derived  for  the  flat  plate,  it  is 
necessary  to  substitute  ^^  *  t^  --  k^  md  r  »  r^  ^  k.    fhese  substltu- 
tlons  are  identicalljr  true  frcsa  the  gecaietry  and  aserely  sei-ve  to  shift 
the  zero  axis  froai  the  tube  centerllne  to  the  tube   v»U.  Making  tte 
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•tibctitutioos  and  simplifying,  we  have 


(68) 


If  tbe  liquid  film  is  thin,  r  »  x  and  r  »  x  .  For  this  condi- 
tioea  considerable  simplification  eiaiy  be  msu^le^  resulting  in  the  follov- 
ing  sfaMir  etrttss  distribution. 

(69) 

This  expression  is  identical  to  £quatic»ii  (2),  which  vas  derived  for 

the  flat  plate,  and  fonos  the  basis  for  the  development  of  the  predicted 

heat  transfer  coefficients. 

In  the  experimental  prc^rasa,  the  Btaxiffiuc  f ilsi  thiclmess  measxxred 
was  of  the  order  of  .02$  inches.     "Riis  represents  %  of  the  tube  radius. 
By  substituting  these  values  in  Equation  (68),  it  is  easily  Bhonax  that 
tlie  maxiimm  error  which  is  intx*oduced  by  the  siiaplifying  asstsiiption 
that  the  film  is  thin  is  c»ily  3*5ic*    TtiXa  maxiniun  error  vill  c^cur  vhen 
considering  the  vail  shear  stress  and  vill  only  be  apparent  i^ien  the 
vapor  shear  stress  forces  are  large  oo^pared  to  the  gz^vity  forces. 
In  this  case,  the  film  thictoesses  are  xexttib.  ssMJUer  than  .O^;^  inches. 

mring  verified  that  tj^^e  flat  plate  equations  tci'  shear  stress 
are  approxiiaately  true  for  the  case  of  a  vertical  tube,  it  is  nov  pos- 
sible to  utilise  the  flat  plate  equations  to  predict  the  film  thick- 

rr  ^ 
«-w-  w.  -  .«^>.-*^..  ^*        *^    .     For  laminar  flov.  Equation  (22)  applies. 

foae  a  giwn  Reynolds  Suraber,  by  ass\xaing  varicais  values  of  x      and 
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solving  for  tv  ,  we  arc  able  to  plot  x   as  a  function  of  *'-'  . 

Similarly,  for  the  case  of  turbulmt  flow  in  tiie  liquid  fila,  for  a 

glToti  Reynolds  Btsab«r,  a  trial  scd  error  doluticaa  of  Equation  (^) 

will  yield  a  value  of  x  " .  Fro©  ti3/e  definition  of  x  "^^  Equation  (56), 

ths  iMLll  shear  stress  can  be  found  in  terms  of  the  film  thijekaess.  x  • 

'     o 

By  equating  this  expression  for  vail  shear  stress  to  that  froa  the 
force  balance.  Equation  (5),  a  formula    for  x„     in  terms  of      t-».  is 
obtained.     This  has  be«n  done  for  botn  ivBrn.yiar  aiod  turbulent  films  at 
the  i^ynolds  Nuabcrs  for  ^diich  f ilOi  thloleness  readings  wtrft  obtained, 
and  forms  the  basis  for  the  theoretical  cunres  sho'un  on  the  plots  of 
measured  data. 

The  proe<^ure  used  In  obtaining  aMsumred  values  of  f  llxa  thi«|BMHW 
is  as  follows  •     Liquid  vns  reclrctilated  vitii  the  puBi;p,  and  the  flow 
rate  adjusted  to  give  the  BeynoLda  Wmimr  desired,    "not  inxap  was  used 
<m  all  of  these  runs,  rather  than  water  directly  from  the  water  isaln 
as  it  was  found  that  fluctuations  in  city  water  pressure  aade  it  impos- 
sible to  obtain  steady  flow  rates.     With  the  wate.^   flow  rate  c^^stant, 
the  throttling  valv<e  on  the  air  sucticwa  line  was  adjusted  to  give  var- 
loos  Isrvels  of  pressure  drc^  throuj^  the  test  section.    At  each  of 
thase  values  of  presstare  drop,  filia  thickness  reading  were  tal&en  wi'Ui 
tba  ttienaster  probe,  followiii^j  the  procedure  described  in  Section  B, 
Part  II.    Befoi-e  each  scries  of  runs,  filia  thlclmese  readings  were 
taJsw  throng  all  three  of  the  probe  ports  in  the  grou^  being  worked 
with  to  instxre  that  filca  thicloiesaes  were  equal  around  the  circimference 
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of  the  txibe  and  that  the  flow  ves  eve»ily  distrilyi^ed .    The  rav  data 
•nd  ccuiputatio£u>  for  iin&Be  nms  axe  lxicl\2ded  in  Appendix  C  aod  are 
shown  plotted  in  Figures  21  through  27. 

The  static  pressure  level  at  points  in  the  test  section  vas  plot* 
ted  ag^alttst  distance,  and  a  straight  line  curve  dra\m  throtjgh  these 
|K>ints.    The  slope  of  this  curve  \«a8  then  taken  as  the  pressure  grad- 
ient in  the  z -direction.    This  pressure  gradient  could  then  be  easily 
coaverted  to  the  vspor  shear  stress,     Tlr-   ,  from  the  formula 
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lampectioo.  of  the  plotted  resnilts  of  the  film  thiclmeas  meaauro- 
mmxta,   Figures  21  to  27,  reveal*  several  interesting  facts.  First, 
tlM  slMipes  of  the   curves  of  x      versus   La^  correspond  closely  with 
the  shapes  of  the  theoretical  cisrves.  In  geoieral,  the  values  of  meas- 
vred  film  thicioiess  lie  belov  the  predicted  valxies,  this  eleviaticsi  be- 
coming more  noticeable  at  higher  liquid  Beynolds  Nvabers*  The  caeaaured 
values  for   l*  »  0  are  eooBistently  higher  than  predicted.  At  hi^ 
valuer  of  i-<    ,   cnperlBMotal  valises  of  x   agree  well  vith  theoretic 
cal  values.  There  si^pcam  to  be  rather  sharp  change  in  curvatui'e  of 
the  Beasuz*ed  val\»a  at  values  of   t^  in  the  nei*i;hborhDQd  of  5  to  10. 

For  U   3  0,  visual  observations  revealed  that  rather  iars©  sur- 
face vavtts  vere  praaent  an.  the  interface  of  the  li<3\iid  layer.  These 
ymw&B  »ate  film  thickn^^s  MNMRsnHumts,  using  the  i&icraaeter  probe  and 
electrical  continuity  technique,  sccwrwhat  uncertain.  The  ohEKieter 
registered  alternate  breaking  and  making  c^  contact  between  these  sur- 
face vsws  and  the  probe.  In  taking  the  readijagB,  an  attempt  \mB  laade 
to  obtain  a  iSfean  vaive  hei^t  raadt&e,  but  this  procedui'e  is  admittedly 
a  crude  one,  &&   surface  tCTUsion  forces  tended  to  maintain  ccmtact  be- 
tween the  probe  tip  and  the  liquid  surface.  Therefore,  it  is  not  sur- 
prising that  these  readings  were  high,  it  being  believed  that  thick- 
nesses, OS  measured,  were  actually  soeae^^mt  nearer  the  crests  of  the 
iwves  than  the  oean  height. 

'  greater  thaa  a3x3Qt  ^,  the  lar^,  slow  i&oving 
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surface  vawB  vcre  replaced  by  mach  ssallGr  smd  more  rapidly  tBOVlog 
irregular itiea .     Indications  on  the  otomoter  ^ve  quite  sharp  cut-off 
points  for  contact  of  tii^  probe  vith  the  surface.     It  is  felt  that  fila 

^.  ^_      _^       ,^  V  are  accurate  to  within 

.0005  inches.     This  is  leased  oaa  the  shaarp  cut-off  points  and  the  coxi- 
•istency  of  the  reSfdizigiB  an.  several  vxexe  under  identical  conditiosiB. 

The  fact  that  meastired  values  of  x      lie  belov  the  predicted  cur- 
ves lea(te  to  the  conclusion  that  the  velocity  distribution  is  actually 
MBrathin^  otter  tiWB  that  aosisiied.     Although  the  theorexicai  turoulent 
curves  lie  above  the  laaiinar  curves,  it  is  believed  that,  for  the  hi^^r 
▼alAies  of     u-  ,  the  fila  is  e»a«atlally  In  turbulent  motion.    This  ccai- 
cluslon  is  based  en  the  appearance  of  the  f  IIqq  and  on  sone  parelicainary 
recruits  obtained  vit^  a  hot  virc  atURGOAter  in  the  llq^uid  film.     The 
flact  that  the  vslocity  distributicax  may  be  something  other  tliari  the 
yniveraal  distribution  of  Prandtl-Hikurmdse  is  not  altogether  surpiis- 
iag.    The  v»e  of  this  velocity  distributicsa  is  certainly  an  over  ext«i- 
aion  of  their  vxark,  -which  ^«s  derived  froa  ntMunir«Bitnts  of  full  pipe 
floT».     aowevex-,  the  actual  values  of  film  thickness  are  at  least  vith- 
in  a  fev  percent  of  the  predicted  results. 

It  is  believed  that  the  change  in  curvature  of  the  aeej&iired  values 
of  film  thlcki^as  is  due  to  Uie  dsoiping  out  of  the  large  surface  waves. 
This  correspcMids  to  visual  observations.     Ihe  general  character  of  the 
flow,  as  observed  through  the  transparent  test  section  nBad  as  folloTiWt 
For     Ur    »  0,  the  flow  -mM  mainly  cLftiar  vith  relatively  larige  aod  reg- 
ular surface  vaines.    lUke  frequency  of  these  mmm  vrnM  low  and  fairly 
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conatant,  being  of  the  order  of  3  cycles  per  seecaid..  As  tlae  air  fiov 
'Mas  increased,  /these  large  surface  WkV9&  vere  daniMd  oat  \mtil,  at 
around  La<  »  5>  tbe  flox^  appesiared  ^ite  clc^ur  and  traxispsrent  •  For 
yet  higher  vjlilues  of  '^v  ^  ver>'  ainall  aurface  i«ves  were  formed,  with 
&  relatively  high  frequency.  It  iu  estisaated  that  the  height  of  these 
waves  vjas  less  than  .CXX)5  incb^s.  It  is  felt  that  thesa  flnall  waives 
Sire  directly  ocamected  vlth  the  onset  of  tyrbulence  in  ths   liquid  filss.. 
It  is  to  be  noted  in  this  respect,  that  these  soall  vaves  could  be 
attained  without  any  air  flow  at  all,  but  required  very  hig^  liquid 
Reynolds  HMnbers.  Hb  wHunireBMszta  of  this  point  were  obtained,  the 
flow  rate  being  beyond  the  capacity  of  the  liquid  flow  laanoffiieter.  It 
is  estimSfted  that  liquid  Eeynolds  Ifumbers  of  the  order  of  aagaitude 
of  6,000  to  8,000  are  required  to  attain  this  condition. 

On  the  basis  of  the  film  thickness  meaiwenients,  no  definite  in- 
foroAticsi  could  be  deduced  as  to  the  onset  of  turbulence,  as  affected 
by  a  rtBQCX  shear  stress.  Ho'Mever,  the  equations  predicting  film  thi^* 
a«M»  pi'o%'»d  to  be  approxJtoately  correct.  It  may  be  concluded  frcm  l^is 
information  that  the  assumed  velocity  dlst£'ibuticai  is  adequate  for  the 
purpose  of  predicting  beat  transfer  eoefficients.  It  may  also  be  con- 
cliscLed  that  the  equations  derived  for  a  flat  plate  vill  be  eq^ially 
8^>plical>le  to  vejrtical  tubes,  as  long  as  the  film  thlclEness  is  azaokll 
ccopared  to  the  tube  radius. 
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la  flurt  1  "Wie  crltcpla  for  datermlniJGg  ttjc  tpansitlori  point  of 
tte  liquid  filjn  fra&  tfanjinar  tx>  turfcruXe^t  flov  ii«re  £iropc«ed.     .iv/bq 
orit«ria,  aXt^boug^  iMUMd  o&  some  expsrim^atal  work,  are  far  tvcxL  being 
vm^UtA,  «qpfteially  in  the  r«gic»  Gt  appreclabls  in^por  siiear  stresses. 
In  faet,  su?  experiffisntsl  work  has  be«n  fouod  Idiat  can  eked.  ligiit  oxx 
this  plMMEMMnaa  (^  tzwuiition,  as  affected  by  a  vi^per  shMur  stress* 
This  section  vill  discuss  tlMft  ppeliAizutrj  v^i^rk  for  deterBtlninf;  the  tram- 
eitiuri  polat  by  aeaas  of  a  hot  wire  SASHiOBMrter .    Ho  qpmtitative  infor- 
atatioci  has  betuk  obtaiiMd,  Inxt  vork  alons  these  lines  is  currently  in 
progriMa  and  iteould  beccxBe  available  in  the  near  future,  reference  (l6). 

The  instrument  is  a  staodaard  hot  vire  ao«MOGMtt«r|  naBufaotured 
by  the  Flov  Corporation  of  CmibTiA0bf  Maaaaehaeetta  •     It  is  beiiig  used 
ia  e<M9raaaor  blade  eaaeade  stuiiee  lii  the  9as  Turbine  Iabc;ratcry  at 
M.  I.  T.    The  ele^n^t  is  a  .000^  inch  diaBMrtesr  tut^sten  wire,  s^^qpported 
by  a  long  shaft.    By  inseartlng  titie  shaft  through  one  of  the  probe  ports 
on  the  test  section,  the  wire  can  be  introduced  into  the  li<suid  film 
on  the  far  side  c^  the  test  section,    ^ith  currents  of  approytoately 
130  milliaKQps,  {good  indicatlosuv  can  be  obtained  on  an  oeeilloacope .  With 
the  relatively  hi^  current,  tih«  i*fire  burns  oaxt  almoet  instantly  if  ex- 
posed to  air.    This  has  been  the  source  of  some  difficulty,  and  it  mmnf 
he  n«e«Mnry  to  modify  the  wire  eVmimt  co*  Ixi^ver  the  currant  to  ijxre- 
vtmt  thla. 

At  l0v  liquid  Reynold*  !f»dt)«rs,  the  flow  mppman  to  be  guite 
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■Mwwytlble  to  i»i»teadin«»«  ommmA  by  vibratlona.     Fr®limtn«ry  runs 
vmrm  mtOm  vhlle  a  X«arge  air  eaa^rmBBear  ^MB  runziijag  In  the  same  build- 
la®,  t^proxim&Uily  100  feat  atmif*    Very  ragular  velocity  fliiotuaticao* 
vara  dbtMrmd  at  lov  Baynolda  WmAteatB,  t«hlch  dlad  out  aa  tba  llq:uid 
flov  «aa  inereaaad.    Tim  rmm  vara  tbesi  rap«atad  at  night,  vhan  the 
coo^xraaaor  i«a  oot  rutmlag^  aad  theaa  fluetuatiocts  \#er«  abaant. 

For  the  case  of  no  absar  atreaa,  laaainar  flow  in  tbe  liquid  film 
axiated  up  to  Reynolda  »*wber  of  6,000  to  3,000.    A  trip  wtra  vaa  in- 
a«rtad  arc»2ad  the  noszle  throat  at  Ute  aatraoce  to  tbe  test  seotion. 
With  tbia  trip  wire,  oocaaioaal  turbulent  fluctuatiooa  vere  obaarvad 
at  th9  ■m*'1'«^y*  flov  rate  tbat  tha  liquid  flow  ■aaoMfter  could  aaaaura. 
It  la  believed  that  atable  lasainar  flow  eclats  v^  to  Heyiotolda  BTanlbara 
of  6,000.    This  obaarvaticai  haa  at  praaent  aot  baaia  verified. 

With  air  flowing  in  the  taat  taction,  a  eona latently  turbulent 
flow  eouM  be  attalaad  with  lower  liquid  Reynolda  Nwbara.     However, 
aa  waa  evident  from  tbe  film  thlekneaa  maasurementa,  the  film  la  ex- 
tresiely  thin  for  tha  regiona  of  apipraciable  vepor  abear  ttreaa,  and 
difficulty  waa  anicoiaatered  In  toeplag  the  hot  wire  caapletely  in  the 
water  film.     It  aay  be  neceaaary  to  uaa  a  mg»9  viscous  fluid  in  order 
to  get  fila  thlckneaaea  of  sufficient  saagaitude  with  which  to  work. 

Qa  tte  baala  af  the  very  alcetcis^  vcark  done  with  this  technique, 
little  of  a  qonatlt&tive  imture  can  be  concluded.     However,  th«  tech- 
nique appaara  to  be  a  powerful  osm^.    Turbulent  indications  ar&  quite 
Barked  on  the  oscilloaecq^  trace,  although  a  lo<wir  bacls^ound  noiae 
Level  would  be  deaireable  to  enable  aora  aaaet  detexT&inatlcm  c£  tha 
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(itlcm  from  ccsopletely  laminar  cooditlone.     It  is  felt  that,  once 
the  minor  problAEBB  hav*  tN»en  solved^  thfs  toot  wire  ioMHKweter  will  yield 
tbxi  aasver  to  Mm  transition  (|uesti(m» 
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In  vjbev  of  the  preaent  lack  c^  imcmlMt^i^  ccxoceimUog  tiMt  b«httvlor 
of  a  falling  liquid  film  frota  a  fluid  dymmiCB  point  of  viev,  it  i« 
recommdML  ttet  the  first  mcgmrUumtal  efforts  be  directed  along  this 
chaimel, 

i^cifieal2y,  it  is  felt  that  it  i»Jould  be  desirable  to  fully  in- 
vest Igftte  the  trszisition  criteria.     'Ibis  should  be  done  vith  at  least 
three  liquids  of  widely  varying  viscosity.     Consideration  siiouid  be 
given  to  varying  surface  tension  also.    TiadQ  recooxaendation  is  based 
OB  tbfS  fact  that,  in  this  type  of  fluid  flow,  these  foroea  aay  ^ite 
possibly  be  the  only  si^ificant  ooes  in  the  radial  dlrectiori.     It  is 
worthy  of  coament  that  a  ourface  tension  ratio  was  the  correlating 
factor  in  the  friction  factors  presented  in  reference  (11). 

So.  i»orking  vith  other  fluids,  it  vould  be  desirable  to  verify 
the  filj»  thioloMMi  —SMTaacnts  presented  in  this  thesis.    By  gsther- 
ing  sufficient  InfiBCMLtian,  equatiosis  itK>re  adeq^tely  describing  tiw 
velocity  distribution  in  the  turbulent  regions  should  be  developed. 

Mnedk  »ore  information  coocering  friction  factors  in  t%fo  phasa 
coc\u'rent  flov  is  neaded.    As  previously  pointed  out,  this  can  be 
easily  obtained  frcm  the  ]^tts«nt  s«t«^  by  addiiag  a  sieter  to  oeasura 
air  flow  rates. 

After  the  ^r*rodyn«Kiics  of  the  liquid  film  have  been  eufficiently 
defined,  the  analysis  should  be  laodif led  to  include  any  nev  or  better 
Infcxnaation  that  ma^y  have  \>9etk  discovered.    Finally  the  analysis  should 
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be  directly  checked  by  investigation  of  an  «9eriiMnt&I  caaimmvr,  in 

yAiich  provisions  should  be  mode  to  detertrs.ine  point  values  end 
values  of  tbe  condensing  side  tieat  transfer  coefficients. 
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Ap-pmVlX  A  "  HUMSBICAL  EVAimTIQN  OT  MALY  3 

Laminar  Portion 

Foir  tbXa  portioai  of  the  coodnuMtte  film,  tha  si^ificaat  Equations 
are  (IB),  (22),  and  (26). 


(i^>) 


(22) 


(26) 


For  glvi»a  values  of  T'>«'  »•««■«  various  values  of  x   lower  thsin  the 

o 

trans  it  icaa  point,  which  ie  given  by  either  Bspa-'tioa  (29)  or  (jS2),    For 
the  aMRiaad  values  of  x     ,  both  z    and    II,    {~fT\  ^'^^  ^  evaluated 

directly.     -—  l-y-J         say  now  be  found  froe  the  given  valtMS  and 

» 
the  coBtputed  values  of  z  . 

Turbulent  Portion 

For  ccurputation  purposes,  E<iuation  (50)  aiay  be  rearranged  into 

tht  BBore  c<»iLvexiient  form 


(5Sa) 
!nie  evaluation  of  this  equation  entails  a  mnerlcal  integration  pro- 
cedure and  is  best  follo>)ed  by  listing  the  steps  involved. 


:)  tew  ,(SS)  %(o^) 


(at) 


(es) 


x%x):Ti:  .\v;4  =7-f^ 


^m  aad^  isvqX  *_x  Y^  iffU'Cjrv  8sfQt«unr  jtuKm   v^    ^  wnim  flvrJqi  nV 


bam  oa»Sm  oav  «wil  JtaMl  94  woa  ^iw         \^-^J  -^r-    .\SJotrUh 
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Step  1  -  For  given  values  of  La>    smd  n^,   asstnm  vwt'lcw* 
of  X      greater  than  that  at  the  transition  point. 

Step  2  -  Evaluate  M  from  Eiiuation  (^^a). 

atap  5  -  Evaluate  x^"  «  k^*  J{x*  r     tJ*  )(i  -  p^p)  .  At  this 

poiat  in  the  (Kxsspat&tixma   It  va»  aeeeeaery  to  aeaiMke  tiiat 

Step  ^  -  Caleulate  F^  frcm  E^paatioo  (^7), 

Stap  5  -  Calculate  F^  froa  fi^uation  (51). 

9tep  6  -  Fom  the  product,  F-Fg  U^j      , 

Step  7  -  Plot  this  pi^oduct  aa  a  fuactioa  of  x  * . 

Step  8  -  At  various  valtwa  of  x^^ ,  tafee  the  area  mder  this  curve 

vlth  a  plaataeter  or  by  othmr  euitable  graphical  neaaa. 

■mie  area  ie  takan  from  the  tranaitioii  pointy  x'  ,  to 


the  ealeeted  values  (^  x  ^.  ^^4 

o  f     *  -      -  ' 

Step  9  -  Caleulat«  («*  -  »  *)  »  -^  - — -r—^^    I      f^fx 

Step  10  -  Add  to  the  neeult  of  step  9  the  distance  to  the  trata» 


iticxa  point,  s^  ,  foiund  in  the  laminar  ealculatiooe  • 

Step  U  -  Calculate  H^^  from  Sfuatioa  (5i*). 

Step  12  -  Hov   iba.  /if")  ^  ;:;   N/Pe. 
If  13/        7* 


* 


Where  z  ia  the  insult  of  step  10. 


Stes>  13  -  Plot  !^  U^j    ^mrvm   »*  aed  M^^^.  There  will  result  a 
double  series  of  cuxrves  with  Tv  aad  IL  as  pararaeterSi 
The  results  of  the  above  calculations  ai*e  ooaatained  in  Table  X  ana  are 
sho-wn  plotted  in  Figures  7  through  12. 
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Theoretical  Predlctlaris  of  Film  Thlcteesg 

l.agii!:Mur     -  Ttm  laaainar  predictions  result  directly  from  Equation 
(22),  which  may  be  rearrar>ged  ixi  the  fcaroi 


* 

For  a  given  K^  ,   the  plot  is  made  by  aseumitig  values  of  x  arid,  solving 
tor      L/v  . 

Turbulent  -  Tne  turbulent  predictions  result  froea  the  followiug 
procedure:  First,  consider  Equation  (^9),  which  relates  rip  a^A  x  , 


Re     o 


For  x_^  >  30 

o  ' 


JL    »  ^Jk   »  (^  .  5  x^^  .  2.5  x/  m  x^^) 
For  5  <  x^'  <  30 

X  -  ^^Jk   .  (12.55  -  8.05  x^^  .  5  x^^  in  x^*) 

For  a  given  TL  ,  this  equation  is  solved  by  trial  and  error  for  x  ^. 
How,  by  eliminating  To  betwee.-i  the  defiaition  of  x  ' ,  >f:quatioQ  (33), 

and  the  force  balaiice.  Equation  (5)»   Tv  may  be  expressed  in  terns 

>      « 
of  X   and  x_  , 

o      o  '  o^ 

With  the  valtae  of  x   determined  ny  the  K^^  ,  and  by  oesiming  various 

values  of  x_  ,  7^   can  be  solved  for. 
o  ' 


iOX- 

.r.-.ci  .:aj  ..i  bsK^iMnman  9tf  x^a  i  (22) 


« 


V 


^,     tot  tor-SB  br&  l&tt^  v.if  f«vr,^«  *)    «.j#«upe   ......   ,   JH  .-wrrls  «  -wt 

« 


mttyftev  SLnlawe&fi  vrf  *^^  ,    .  ^-i»^f^i>      j^  ^^  wrx«v  ««is; 


o 


RJiK 


o 
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Ti»  above  procediires  for  both  the  laminar  and  turbulent  films  will 
yield  the  predicted  filjs  thiekneases  sho^m  plotted  in  Figures  21  to  27. 
Measured  V&luea  of  Film  Thiekmess 

Measured  data  is  contained  in  I^le  m.  Presetire  readings  were 
plotted  aad  the  straight  line  slopes  taken  to  give  dp/dz.  Froci  Equa- 
tic«  (71),  by  neglecting  p  ccrapared  to  p,  Tv  caa  r«adily  be  coirputed. 
Sl&ilarly  x   can  be  found  fraa  its  definition  aad  the  measured  values 
of  X  .  Calculaticsis  are  aJjto  included  in  l^le  lil. 


-SOI- 


JUlv  amixi 

•  TS  o»  IS  ■»'Xi/8iI  aX  b^ 

aawxav  harxmjum  c^tf  du«  csol 


:ite  «»a«» 


•IT 


M*^ 


ftitmMfntrfT    ^itF?'?   '>/•■>    t.*..-rr«  '    ^*»../e^^uJ 


o 


-103- 
^  III 


LiauiD  n     m  538 


F.e 


Vteter  Orifice:     0.10  lii.  Hg. 

am        .095 

Uater  Taap:     75.5°  F. 


Wrter  Vl»co6lty:     1.955  X  lo"^  Ib-Mc/ft^ 

"f. 

55B 

•■ 

Of 

Probe 

568    l/ln. 
vilall  Heading: 

.5265 

iachsa  HgO 

below  atta. 
P3          P4 

P5 

»6 

in.  HgO 
ft 

dp/dz       tT 

Film 

Probe 

Rding 

itxbtam 
0            0 

0 

0 

0 

0 

0 

0 

0 

0 

.502 

24.5 

15.90 

0 

0 

.16 

.22 

.24 

.55 

.090 

1.06 

.5095 

17.0 

9.65 

0 

0 

.20 

.35 

.42 

.54 

.120 

1.42 

.511 

15.5 

8.80 

0 

0 

.41 

.42 

.80 

1.7 

.220 

2.60 

.5165 

uo.o 

5.68 

0 

.07 

.55 

1.5 

1.5 

2.0 

.413 

4.94 

.5178 

3.7 

4.94 

.05 

M 

1.0 

1.85 

2.15 

2. 85 

.588 

6.95 

.5192 

7.3 

4.15 

.^5 

i.k 

2.9 

3.7 

4.4 

5.15 

1.175 

13.90 

.5209 

5.6 

3.18 

1.5 

2.7 

4.5 

5.9 

7.5 

8.2 

1.575 

18.65 

.5216 

4.9 

2.78 

2.7 

>.25 

6.9 

a.75 

10.5 

11.7 

2.175 

25.70 

.5221 

4.5 

2.55 

4.5 

6.6 

10.2 

15.1 

15.4 

17.1 

5.150 

57.30 

.5225 

4.2 

2.58 

6.9 

9.9 

14.9 

18.5 

21.6 

24.5 

4.200 

49.60 

.5229 

5.6 

2.04 

U.55 

16.1 

25.55 

29.2 

54.3 

58.9 

6.550 

74.90 

.5251 

3.4 

1.93 

-eoi- 


^x,  ^pa 


OUk* 


^•«-» 


tc- 


•flrfOOl 


out 


5d»^      0»TI    {^i^-     ^''^•i       ji^v.       c 


-■•Ui 


.-     o  r 


>i.<" 


■,4 

1 

CC.i. 

^.4 

t 

<k),- 

e?.i 

i-'^. 

^o 


l.j^J.         vl. 


-<*-'■ 


.ni\x     tvic         <ifH>f) 


>flB^4l  lfOd«Cf  OgH 


4« 


^4.  <^^.  y.rv. 

T.X      06.        S4.        X4, 


^^'C 


:^.Gi. 


s« 


0 
TO. 

'.£        4.i 

r.s 


0 
0 
0 
0 
0 

to. 

.XI 


Liquid  K 


He 


k^^ 


-104- 
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